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A  US!  RACT 

New  pii ra mete ri /a lions  for  ihe  spectral  dissipation  of  wind  generated  waves  are  proposed.  ITie  rates  of 
dissipation  have  no  predetermined  sj>ectral  shapes  and  are  functions  of  the  wave  spectrum  and  wind  speed 
and  direction,  in  a  way  consistent  with  observations  of  wave  breaking  and  swell  dissipation  properties. 
Namely,  the  swell  dissipation  is  nonlinear  and  proportional  to  the  swell  steepness,  and  dissipation  due  to  wave 
breaking  is  nonzero  only  when  a  nondimensional  spectrum  exceeds  the  threshold  at  which  waves  sire  observed 
to  start  breaking.  An  additional  source  of  short wave  dissipation  is  introduced  to  represent  the  dissipation  of 
short  waves  due  to  longer  breaking  waves.  A  reduction  of  the  wind  wave  generation  of  short  waves  is  meant  to 
account  for  the  momentum  flux  absorbed  by  longer  waves.  These  paramcicrizations  arc  combined  and 
calibrated  with  the  discrete  interaction  approximation  for  the  nonlinear  interactions.  Parameters  arc  adjusted 
to  reproduce  observed  shapes  of  directional  wave  spectra,  and  the  variability  of  spectral  moments  with  wind 
speed  and  wave  height  I  lie  wave  energy  balance  is  verified  in  a  wide  range  of  conditions  and  scales,  from  the 
global  ocean  to  coastal  settings.  Wave  height,  peak  and  mean  periods,  and  spectral  data  are  validated  using  in 
situ  and  remote  sensing  data.  Some  systematic  defects  arc  still  present,  but.  overall,  the  paramcicrizations 
probably  yield  the  most  accurate  estimates  of  wave  parameters  to  date.  Perspectives  for  further  improvement 
are  also  given, 


1.  Introduction 

a.  On  phase-averaged  models 

Spectral  wave  modeling  has  been  performed  for  the 
last  50  years,  using  the  wave  energy  balance  equation 
(Gelci  el  al.  1957),  This  approach  is  based  on  a  spectral 
decomposition  of  the  surface  elevation  variance  across 
wavenumbers  A  and  directions  0,  The  spectra  density  F 
evolves  tn  live  dimensions  that  arc  the  two  spectral 
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dimensions  k  and  0.  the  two  physical  dimensions  of  the 
ocean  surface  (usually  longitude  and  latitude),  and  time  i: 


where  the  Lagrangian  derivative  is  the  rale  of  change  of 
the  spectral  density  following  a  wave  packet  al  tls  group 
speed  in  both  physical  anil  spectral  spaces.  In  particular, 
this  spectral  adveclion  includes  changes  in  direction  due 
to  the  earth's  sphericity,  as  well  as  refraction  over  varying 
topography  (e.g..  Munk  and  Traylor  1947;  Magne  et  al. 
2(X)7)  and  currents,  and  changes  in  wavelength  or  period 
in  similar  conditions  (Barber  1949). 

The  spectral  source  functions  on  the  right-hand  side  of 
Eq.  (1)  are  grouped  into  iheir  atmospheric  (.VlIm).  non¬ 
linear  scattering  (Sn i).  oceanic  (,Sf>c).  and  bottom  (.S^,) 
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sources.  This  grouping,  like  any  other,  is  largely  arbitrary. 
For  example,  waves  that  break  are  nonlinear.  Thus,  the 
effects  of  breaking  waves,  which  are  contained  in  S1K, 
is  intrinsically  related  to  the  nonlinear  evolution  term 
contained  in  Snl.  Yet,  compared  to  the  usual  separation 
of  deep-water  evolution  by  wind  input,  nonlinear  inter¬ 
actions,  and  dissipation  (e.g..  WAMDI  Group  1988),  this 
grouping  has  the  benefit  of  identifying  where  the  energy 
and  momentum  are  going  to  or  coming  from,  which  is 
a  necessary  feature  when  ocean  waves  are  used  to  drive 
or  are  coupled  with  atmospheric  or  oceanic  circulation 
models  (e.g..  Janssen  el  al.  2004;  Ardhuin  el  al.  2008b). 

Here,  5aIfn,  which  gives  the  llux  of  energy  from  the 
atmospheric  nonwave  motion  to  the  wave  motion,  is  the 
sum  of  a  wave  generation  term  Sm  and  a  wind-generation 
term  Sollt  (often  referred  to  as  negative  wind  input,  i.e., 
a  wind  output).  The  nonlinear  scattering  term  Sn\  repre¬ 
sents  all  processes  that  lead  to  an  exchange  of  wave  en 
ergy  between  the  different  spectral  components.  In  deep 
and  intermediate  water  depth,  this  is  dominated  bv  cubic 
interactions  between  quadruplets  of  wave  trains,  while 
quadratic  nonlincarilies  play  an  important  role  in  shal¬ 
low  water  (e.g.,  WISE  Group  2007).  The  ocean  source 
SiK  may  accommodate  wave-current  interactions 1  and 
interactions  of  surface  and  internal  waves.  The  oceanic 
source  term  5<K  is  restricted  to  wave  breaking  and  wave- 
turbulence  interactions. 

The  basic  principle  underlying  Eq.  (1)  is  that  waves 
essentially  propagate  as  a  superposition  of  almost  linear 
wave  groups  that  evolve  on  longer  lime  scales  as  a  result 
of  weak-in-the-mean  processes  (e.g.,  Komen  el  al  1994). 
Recent  reviews  have  questioned  the  possibility  of  fur 
liter  improving  numerical  wave  models  without  chang¬ 
ing  this  basic  principle  (Cavaleri  2006).  Although  this 
may  be  true  in  the  long  term,  we  demonstrate  here  that  it 
is  possible  to  improve  our  model  results  significantly  by 
including  more  physical  features  in  the  source  term  pa- 
ramelenzalions.  The  main  advance  that  we  propose  is  the 
adjustment  of  a  dissipation  function  without  any  pre¬ 
scribed  spectral  shape,  based  on  our  empirical  knowledge 
of  the  breaking  of  random  waves  (Banner  cl  al.  2000; 
Babanin  el  al.  2001)  and  the  dissipation  of  swells  over 
long  distances  (Ardhuin  el  al.  2009a).  The  present  for¬ 
mulations  are  not  based  on  a  detailed  physical  model  of 
dissipation  processes,  but  they  demonstrate  that  progress 
is  possible.  This  effort  opens  the  way  for  physical  pa- 
ramelerizations  (e.g..  Filipoi  el  al  2010)  that  will  even¬ 
tually  provide  new  applications  for  wave  models,  such  as 


1  In  I  lie  presence  of  variable  current,  I  lie  source  of  energy  for  ihe 
wave  field  is  Ihe  work  of  [lie  radial  ion  stresses.  Ii  is  nol  explicit 
when  Ihe  energy  balance  is  written  as  an  action  balance  (eg,  Komen 
ct  al  1994). 


the  estimation  of  statistical  parameters  for  breaking 
waves,  including  whilecap  coverage  and  foam  thick¬ 
ness.  Other  efforts,  though  less  empirical  in  nature,  are 
also  under  way  to  arrive  al  belter  paramelcrizalions  (e.g.. 
Banner  and  Morison  2010;  Babanin  et  al.  2(X)7;  Tsagaieli 
2008),  but  they  have  yet  to  produce  a  practical  alternative 
for  wave  forecasting  and  hindcasting. 

b.  Shortcomings  of  existing  pa  ram  e  teri  z  n  tions 

Until  the  work  of  van  dor  Westhuysen  el  al.  (2(X)7),  none 
of  the  wave  dissipation  paramelerizalions  presented  in 
the  literature  had  a  quantitative  relationship  with  observed 
features  of  wave  dissipation  but  rather  they  were  adjusted 
to  close  the  wave  energy  balance.  Komen  el  al.  (1984. 
hereinafter  KHH)  have  produced  a  family  of  paramc- 
terizations  loosely  justified  by  the  "random  pulse”  theory 
of  Hasselmann  (1974).  In  deep  water,  these  parameter! 
zations  take  a  generic  form: 


0)  -  CdJ"krHt 


k 


/'(*.  0 ). 
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in  which  Ct\s  is  a  negative  constant,  F(A,  0)  is  the  spectral 
density,  and  If  —  4  v/Xf  M)dkdQ  is  the  signilicant 
wave  height.  The  energy  weighted  mean  wavenumber 
kr  is  defined  from  the  entire  spectrum  as  follows: 
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where  r  is  a  real  constant.  Usual  choices  are  r  —  -0.5 
(WAMDI  Group  1988)  or  r  ~  0.5  (KHH:  Bidlol  el  al. 
2005). 

These  parameterizalions  are  still  widely  used  in  spite 
of  inconsistencies  in  the  underlying  theory  Indeed,  if 
whitccaps  do  act  as  random  pressure  pulses,  their  aver¬ 
age  work  on  the  underlying  waves  only  occurs  because 
of  a  phase  correlation  between  the  vertical  orbital  ve¬ 
locity  Held  and  the  moving  whilecap  position,  which 
travels  with  the  breaking  wave.  In  reality  the  horizontal 
shear  is  likely  to  be  the  dominant  mechanism  (Longuel- 
Higgins  and  Turner  1974),  but  the  question  of  correla¬ 
tion  remains  the  same.  For  any  given  whilecap.  such  a 
correlation  cannot  exist  for  all  spectral  wave  compo¬ 
nents:  a  whiteeap  that  travels  with  one  wave  leads  to  the 
dissipation  of  spectral  wave  components  that  propagate 
in  similar  directions,  with  comparable  phase  velocities. 
However,  whilecaps  moving  in  one  direction  will  give 
(on  average)  a  zero  correlation  for  waves  propagating  in 
the  opposite  direction  because  the  position  of  the  crests 
of  these  opposing  waves  are  completely  random  with 
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respect  to  the  whilecap  position.  As  a  result,  not  all  wave 
components  are  dissipated  hv  a  given  whitccap  (others 
should  even  be  generated),  and  the  dissipation  function 
cannot  take  the  spectral  form  (2). 

A  strict  interpretation  of  the  pressure  pulse  model  gives 
a  zero  dissipation  for  swells  in  the  open  ocean  because  the 
swell  wave  phases  are  uncorrelated  to  those  of  the  shorter 
breaking  waves.  There  is  only  a  negligible  dissipation  due 
to  short-wave  modulations  by  swells  and  preferential 
breaking  on  the  swell  crests  (Phillips  1963;  Hasselmann 
1971:  Ardhuin  and  Jenkins  2005).  Still,  the  dissipation 
functions  given  by  Eq.  (2)  are  applied  to  the  entire  spec¬ 
trum,  including  swells,  w  ithout  any  physical  justilicalion. 

In  spile  of  its  successful  use  for  the  estimation  of  the 
significant  wave  height  //v  and  peak  period  Tr ,  these  fixed- 
shape  parameteri/ations,  from  KHH  up  to  Bidlot  el  al. 
(2007a),  have  built-in  defects  that  limit  the  accuracy  of 
the  mode  results.  Most  conspicuous  is  the  spurious  am¬ 
plification  of  wind  sea  growth  in  the  presence  of  swell 
(c.g.,  van  Vleddcr  and  Hurdle  2(X)2),  which  is  contrary  to 
all  of  the  observations  (Dobson  ct  al.  1989;  Violantc- 
Carvalho  et  al.  2(X)4;  Ardhuin  el  al.  2007). 

Also,  these  parameteri/ations  typically  give  a  decreasing 
dissipation  of  swell  with  increasing  swell  steepness,  con¬ 
trary  to  all  of  the  observations  from  Darbyshire  (1958)  to 
Ardhuin  el  al.  (2(X)9b).  1  his  effect  is  easily  revealed  by 
taking  a  sea  stale  composed  of  a  swell  of  energy  E\  and 
mean  wavenumber  k\  and  a  windsca  of  energy  E2  and 
mean  wavenumber  k2.  with  k2  >  k\.  The  overall  mean 
wavenumber  is* 

kr2E2)l(E\  +  E2)]l  r.  (4) 

Al  low  frequency  (i.e.,  small  wavenumbers)  the  dissi¬ 
pation  lor  a  given  value  of  k  is  given  by  the  lirst  term, 
proportional  to  k,  5(/T|  +  E2).  Now,  if  we  keep  k\.  k2  an 
E2  constant  and  only  increase  the  swell  energy  E\*  the 
relative  change  in  dissipation  is,  according  to  (2),  pro¬ 
portional  to  v  =  3.5[ (A ! /Av)r  -  1  ]/r  +  2.  For  r  =  0.5,  as 
used  by  Bidlot  et  al.  (2005,  hereinafter  BAJ),  v  is  neg¬ 
ative  (i.e.,  the  dissipation  decreases  w  ith  increasing  swell 
energy)  for  k\/kr  <  0.51.  For  equal  energy  in  sea  and 
swell,  this  occurs  when  k\/k2  <0.3,  which  is  generally  the 
case  with  sea  and  swell  in  the  ocean.  This  erroneous  de¬ 
crease  of  swell  dissipation  w  ith  increasing  swell  steepness 
is  reduced  when  the  model  frequency  range  is  limited  to 
a  maximum  frequency  of  0.4  Hz,  in  w  hich  case  the  lowest 
winds  (less  than  5ms  )  are  unable  to  produce  a  realistic 
wind  sea  level,  hence  limiting  the  value  of  k,  to  relatively 
small  values. 

An  alternative  and  widely  used  formulation  has  been 
proposed  by  Tolman  and  Chalikov  (1996),  and  some  of 


its  features  are  worth  noting.  We  refer  to  this  param¬ 
eterization  as  TC,  including  the  later  adjustments  by 
Tolman  (2002b).  It  combines  two  distinct  dissipation 
formulations  for  high  and  low  frequencies,  with  a  tran¬ 
sition  al  2  times  the  wind  sea  peak  frequency.  Janssen 
ct  al.  (1994)  had  already  introduced  the  use  of  two  terms. 
k  and  k “  in  Eq  (2),  to  represent  different  balances 
in  the  high  and  low-frequency  parts  of  the  spectrum. 
However,  both  terms  are  still  multiplied  by  the  same 
factor,  C^sg0  T?  V/s4<5.  In  TC,  these  two  dissipation 
terms  arc  completely  distinct,  the  low-frequcncv  part 
being  a  linear  function  of  the  spectral  density  and  pro¬ 
portional  to  wind  friction  velocity,  //*:  the  high-frequency 
part  is*  also  linear  and  proportional  to  //;.  In  TC  the 
frequency  dependence  of  the  two  terms  is  also  pre¬ 
scribed.  Tolman  and  Chalikov  (1996)  further  included 
swell  attenuation  bv  the  wind,  here  noted  Soul,  based 
on  numerical  simulations  of  the  airflow  above  waves 
(Chalikov  and  Belevich  1993).  At  relatively  short 
fetches,  the  5in  and  SiK  terms  in  TC  are  typically  a  factor 
of  2-3  limes  smaller  than  those  of  Janssen  et  al.  (1994). 
This  small  magnitude  of  the  source  terms  was  found  to 
produce  important  biases  in  wave  growth  and  wave 
directions  at  short  fetch  (Ardhuin  et  al.  2007).  Another 
successful  set  of  parameteri/ations,  for  high  winds  con¬ 
ditions,  is  the  combination  by  Makin  and  St  am  (2003),  but 
it  does  not  produce  accurate  results  in  moderate  sea  slates 
(Lcfcvrc  cl  al.  2004).  Polnikov  and  Inocenti.nl  (2(X)S)  have 
also  proposed  new  source  term  formulations,  but  their 
results  appear  to  generally  be  less  accurate  than  those 
with  the  model  presented  here,  in  particular  for  mean 
periods. 

Considering  the  observed  strong  wave  height  gradi¬ 
ents  in  rapidly  varying  currents,  Phillips  (1984)  proposed 
a  dissipation  rate  proportional  to  the  nondimensional 
spectrum,  also  termed  the  saturation  spectrum  Banner 
et  al.  (2002)  indeed  found  n  correlation  of  the  direction- 
integrated  saturation  B  with  the  breaking  probability  of 
dominant  waves.  In  particular,  they  found  that  breaking 
occurs  when  B  exceeds  a  threshold  Br.  Alves  and  Banner 
(2(X)3)  proposed  to  deline  the  dissipation  by  BIBr 
to  some  power,  multiplied  by  a  Konicn-type  dissipation 
term.  Although  this  approach  avoided  the  investigation 
of  the  dissipation  of  nonbreaking  waves,  it  imported  all 
of  the  above-mentioned  defects  of  the  KHH  parame- 
lerizations.  Further,  theii  value  for  /?r,  which  was  much 
higher  than  suggested  by  the  observations,  tends  to 
disconnect  the  parameterization  from  the  observed  ef¬ 
fects  (Babanin  and  van  der  Wcslhuysen  2(X)8). 

The  use  of  a  saturation  parameter  was  taken  up  again 
by  van  der  Wcslhuysen  et  al.  (2(X)7,  hereafter  WZB). 
They  followed  Banner  ct  al.  (2(X)2)  and  integrated  the 
saturation  spectrum  over  directions: 
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5(A)  =  kT(k.  0')<I0’.  (5) 
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From  this,  the  WZB  source  function  was  delined  as 


^  ^lUAV V.H  v 


T(A\  tf),  (6) 


basins  larger  than  1 000  km.  Because  of  the  differences 
in  coastal  and  larger-scale  sea  stales  (e.g.,  Long  and 
Resio  2007),  verifying  the  source  function  parameteri¬ 
zation*  at  all  scales  is  paramount,  in  order  to  provide 
a  robust  and  comprehensive  parameterization  of  the 
wave  dissipation. 

c.  A  new  set  of  parameter!  zations 


where  Ct|*.w7.n  >s  a  negative  constant,  Br  is  a  constant 
saturation  threshold,  and  p  is  a  coefficient  that  varies 
both  with  the  wind  friction  velocity  //*  and  the  degree 
of  saturation  B{k)IBr  with,  in  particular,  p  **  0  for 
B(k )  <  0.8/?,..  For  nonbreaking  waves,  when p  **  0,  the 
dissipation  is  too  large  by  at  least  one  order  of  magni¬ 
tude.  making  the  parameterization  unlit  for  occanic- 
scale  applications,  with  wave  heights  in  the  Atlantic 
underpredictcd  by  about  50%  (Ardhuin  and  Le  Boyer 
2006),  Van  dcr  Westhuvsen  (2007)  eventually  replaced 
this  dissipation  of  nonbreaking  waves  with  the  Komen- 
lypc  form  proposed  by  the  WAMD1  Group  (1988). 

Further,  the  increase  in  p  with  the  inverse  wave  age, 
uJC.  tends  to  increase  S (K  at  high  frequency,  and  was 
needed  to  obtain  a  balance  with  the  Satm  term  in  Eq.  ( 1 ). 
This  indicates  that,  in  addition  to  the  value  of  the  satu¬ 
ration  /?r,  other  factors  may  be  important,  such  as  the 
directionality  of  the  waves  (Banner  el  al.  2002).  Other 
observations  clearlv  show  that  the  breaking  rate  of  high- 
frequency  waves  is  much  higher  for  a  given  value  of  /?, 
probably  due  to  the  cumulative  effects  bv  which  the 
longer  waves  arc  modifying  the  dissipation  of  shorter 
waves. 

Banner  ct  al.  (1989)  and  Melville  el  al.  (2002)  have 
shown  that  breaking  waves  suppress  the  short  waves  on 
the  surface.  Young  and  Babanin  (2006)  arrived  al  the 
same  conclusion  from  the  examination  of  pre-  and  post- 
breaking  wave  spectra,  and  proposed  a  parameterization 
for  StK:  that  included  a  new  '‘cumulative  term"  to  repre¬ 
sent  this  effect.  Yet.  their  estimate  was  derived  for  very 
strong  wind-forcing  conditions  only.  Their  interpretation 
of  the  differences  in  parts  of  a  wave  record  w  ith  break¬ 
ing  and  nonbreaking  waves  implies  an  underestimation 
of  the  dissipation  rales  because  the  breaking  waves  have 
already  lost  some  energy  when  they  are  observed  and 
the  nonbreaking  waves  arc  not  going  to  break  right  after 
they  have  been  observed.  Also,  since  the  spectra  are 
different,  nonlinear  interactions  must  be  different,  even 
on  this  relatively  small  time  scale  (e.g..  Young  and  van 
Vledder  1993,  Fig.  5).  and  the  differences  in  spectra  may 
not  he  the  result  of  dissipation  alone 

Finally,  the  recent  measurement  of  swell  dissipation 
by  Ardhuin  cl  al.  (2009a)  has  revealed  that  the  dissipa¬ 
tion  of  nonbreaking  waves  is  essentially  a  function  of  the 
wave  steepness,  and  a  very  important  process  for  ocean 


It  is  thus  lime  to  combine  the  existing  knowledge  on 
the  dissipation  of  breaking  and  nonbreaking  waves  to 
assemble  a  parameterization  for  the  dissipation  of 
waves.  Our  objective  is  to  provide  a  robust  parameter¬ 
ization  that  improves  on  existing  wave  models.  For  this 
we  will  use  the  parameterization  bv  B  AJ  as  a  benchmark 
because  it  was  shown  to  provide  the  best  forecasts  on 
global  scales  (Bidlot  el  al.  2007b)  before  the  advent  of 
the  parameterization*  presented  here.  BAJ  is  also  fairly 
close  to  the  widely  used  Wave  Action  Model  (WAM) 
“cycle  4”  parameterization  employed  by  Janssen  el  al. 
(1994). 

We  will  First  present  a  general  form  of  the  dissipation 
terms  based  on  observed  wave  dissipation  features.  The 
degrees  of  freedom  in  the  parameterization  are  then 
used  to  adjust  the  model  results.  A  comprehensive  val 
idation  of  wave  parameters  is  presented  using  field  ex¬ 
periments  and  a  1-yr  hindcast  of  waves  al  the  global  and 
regional  scales,  in  which  all  existing  wave  measurements 
are  considered,  with  significant  wave  heights  ranging 
from  0  to  17  m.  The  model  is  further  validated  with  in 
dependent  data  at  regional  and  global  scales.  Model 
hindcasls  and  forecasts  al  global  and  regional  scales  are 
available  online  (http  //tiny  url.com/yetsofy),  covering  at 
least  the  years  2002-10. 

Tests  and  verification  in  the  presence  of  currents,  and 
using  a  more  realistic  parameterization  of  wave-wave 
interactions  will  be  presented  in  future  publications. 
These  mav  also  include  some  replacement  of  the  arbi¬ 
trary  choices  made  here. 

2.  Parameterization* 

In  this  section  we  present  the  results  of  the  integration 
of  the  energy  balance.  Because  numerical  choices  can 
have  important  effects  (e.g.,  Tolman  1992;  Hargreaves 
and  Annan  2000),  a  few  details  should  be  given.  All 
calculations  are  performed  within  the  WAVEWATCH 
111  modeling  framework  (Tolman  2008,  2009;  hereinaf¬ 
ter  WWATCH),  using  the  third-order  spatial  and  spec¬ 
tral  adveclion  scheme  and  including  modifications  of 
the  source  terms  described  here  The  source  terms  are 
integrated  with  the  fully  implicii  scheme  of  Hargreaves 
and  Annan  (2000),  combined  with  the  adaptation  lime 
step  and  limiter  method  of  Tolman  (2(X)2a),  in  which  a 
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minimum  lime  step  of  10  s  is  used,  so  that  the  limiter  on 
wind-wave  growth  is  almost  never  activated.  The  di 
agnostic  tail,  proportional  to  /  \  is  only  imposed  at  a 
cutoff  frequency  /  set  to 

4  4  m  4»r  ^ 

Here,  we  lake  /F m  =  10  and  define  the  mean  frequency 
os  f„  =  l/rm(i  Hence.  /-  is  generally  above  the  maxi¬ 
mum  model  frequency  that  we  fixed  at  0.72  Hz,  and  the 
high-frequency  tail  is  left  to  evolve  freely.  Some  com¬ 
parison  tests  are  also  performed  with  other  parameter¬ 
ization  using  a  lower  value  of  /.,  typically  set  at  2.5/, , 
(Bidlot  ct  al.  2007a).  In  such  calculations,  although  the 
net  source  term  may  be  nonzero  at  frequencies  above / , 
there  is  no  spectral  evolution  due  to  the  imposed  tail. 

a.  Nonlinear  wave-wave  interactions 

All  the  results  discussed  and  presented  in  this  section 
are  obtained  with  the  discrete  interaction  approxima¬ 
tion  (DIA)  of  Hassclmann  et  al.  (1985).  The  coupling 
coefficient  that  gives  the  magnitude  of  the  intcraetions 
is  CnI.  Based  on  comparisons  with  exact  calculations. 
KHH  adjusted  the  value  of  Cni  to  2.78  X  10  ,  which  is  the 
value  used  by  Bidlot  et  al.  (2005).  Here,  this  constant 
will  be  modified  slightly.  The  DIA  parameterization  is 
well  known  for  its  shortcomings  (Banner  and  Young 
1994),  and  the  adjustment  of  other  parameters  probably 
compensates  for  some  of  these  errors.  This  matter  will  be 
fully  discussed  in  a  future  publication. 

b.  Swell  dissipation 

Observations  of  swell  dissipation  are  consistent  with 
the  effects  of  friction  at  the  air-sea  interface  (Ardhum 
et  al.  2009a),  resulting  in  a  11  ux  of  momentum  from  the 
wave  held  to  the  wind  (Harris  1966).  We  thus  write  the 
swell  dissipation  as  a  negative  contribution  Soul,  which  is 
added  to  S,n  to  make  the  wind- wave  source  term  Salm. 

Using  the  method  of  Col  lard  el  al.  (2009).  a  systematic 
analysis  of  swell  observations  by  Ardhum  et  al.  (2009a) 
showed  that  the  swell  dissipation  is  nonlinear,  which  is 
possibly  related  to  a  laminar-to-turbulent  transition  of 
the  oscillatory  boundary  layer  over  swells.  Defining  the 
boundary  Reynolds  number  Re  -  4iflirh!/orb/*'</«  where 
//(>rh  and  rtn,h  are  the  significant  surface  orbital  velocity 
and  displacement  amplitudes,  and  vl}  is  the  air  viscosity, 
we  take,  for  Re  less  than  a  critical  value  Re,.. 

S„Jk.  0)  =  -cJ±H2ks/2v<T)F(k,  0).  (8) 

where  the  constant  Ctjs%  is  equal  to  1  in  Dore  (I978)s 
laminar  theory  and  the  radian  frequency  rr  is  related  to  k 
via  the  dispersion  relation. 


When  the  boundary  layer  is  expected  to  be  turbulent, 
for  Re  >  Re,.,  we  take 

SM(k.»)  =  -^(WyUoit,/g)F(k.^.  (l>) 

‘  H 

A  few  tests  have  indicated  that  a  threshold  Re,.  =  (2  X 
10  m )///,  provides  reasonable  results,  although  it  may 
also  be  a  function  of  the  wind  speed,  and  we  have  no  ex¬ 
planation  for  the  dependence  on  II,.  A  constant  threshold 
close  to  2  X  105  provides  similar — but  less  accurate — 
results.  Here,  we  shall  use  CJsv  ”  I  2. 

The  parameterization  of  the  turbulent  boundary  layer 
is  a  bit  more  problematic  and,  in  the  absence  of  direct 
measurements  in  the  boundary  layer,  leaves  room  for 
speculations.  Assuming  a  constant  /.  in  (9),  Ardhum 
et  al,  (2009b)  found  that  swell  observations  are  consis¬ 
tent  with  0.004  <  /.  <  0.013.  From  the  analogy  with  an 
oscillatory  boundary  layer  over  a  fixed  bottom  (Jensen 
et  al.  1989),  these  values  correspond  to  a  surface  with 
a  very  small  roughness.  Because  we  also  expect  the  wind 
to  influence /.,  we  include  the  adjustable  effects  of  wind 
speed  on  the  roughness,  and  an  explicit  correction  ol  /. 
in  the  form  of  a  Taylor  expansion  to  first  order  in  uju,nw 

fr=Sl{f,Xi  M  +  H5.»l  +  COS(W  -  0J\~  -).  (10) 

f  "oi+J 

where  /  um  is  the  friction  factor  given  by  Grant  and 
Madsen's  (1979)  theory  for  rough  oscillatory  boundary 
layers  without  a  mean  flow.  Adequate  swell  dissipation 
is  obtained  with  constant  values  of  /,  in  the  range  0.004— 
0.007,  hut  these  do  not  necessarily  produce  the  best  re¬ 
sults  when  comparing  wave  heights  to  observations.  Based 
on  the  simple  idea  that  most  of  the  air-sea  momentum 
flux  is  supported  by  the  pressure  slope  correlations  that 
give  rise  to  the  wave  field  (Donelan  1998;  Pcirson  and 
Banner  2003),  we  have  set  the  surface  roughness  seen  by 
the  oscillatory  flow  c<)  to  .»  small  fraction  of  that  seen  by 
the  mean  flow 


Here.  rA)  is  here  set  to  0.04.  As  a  result,  typical  values  of 
aor\Jz[)  are  as  large  as  2  X  1(I\  with  of  the  order  of 
0.003. 

The  coefficients  xi  and  s\  of  the  0(ujnar i,)  correction 
have  been  adjusted  to  0.018  and  0.015,  respectively, 
the  former  negative  value  giving  a  stronger  dissipation 
for  swells  opposed  to  winds,  when  cos(tt  0U)  <  0.  This 
gives  a  range  of  values  of  /,  consistent  with  the  obser¬ 
vations.  and  reasonable  hindcastsof  swell  decay  (Fig.  1 ), 
with  a  small  underestimation  of  dissipation  for  sleep 


1922 


JOUKNAl  OF  PHYSICAL  OCEANOdRAPhh 


Voi.KMl  40 


3.5 


2.5 


2  15 

X 


0.5 


+ 

+ 


+ 

* 


V 


+ 

♦ 


% 


+ 

+ 

W 


5000 

3.5 
3r 

2.5  * 


6000  7000  8000  9000  10000 

Distance  along  great  circle  (km) 


+  SAR  ENVISAT 
V  WW3  TEST437 
■  ■  No  dissipation 


%  1.5 

</5 

2  i 

X 

0.5 


*  * 


—  + 


2000  4000  6000  8000  10000  12000  14000 

Distance  along  great  circle  (km) 


I  Ki.  1  Comparison  of  modeled  swell  significant  heights,  fol¬ 
lowing  tliL  propagation  of  the  two  swells  shown  by  Ardhuin  el  al. 
(2009a)  with  peak  periods  of  15  s  and  (top)  high  or  (bottom)  low- 
dissipation  rates,  known  biases  in  the  level  2  data  have  been  cor¬ 
rected  following  Collard  et  al.  (2009). 


swells.  Because  an  increase  of  s  i  from  0.8  to  1 . 1  produces 
negative  biases  on  II,  of  the  order  of  30%  at  all  oceanic 
buoys,  the  magnitude  of  the  swell  dissipation  cannot  be 
much  larger  is  than  chosen  here. 

Further  validation  of  the  swell  dissipation  is  provided 
by  the  global-scale  liindcasts  in  section  4. 

c.  Wave  breaking 

Observations  show  that  waves  break  when  the  orbital 
velocity  at  their  crest.  Ut.  comes  close  to  the  phase 
speed,  C.  with  a  ratio  II JC  >  0.8  for  random  waves 
(Tulin  and  Landrim  2001;  Stansell  and  MacFarlanc 
2002;  Wu  and  Ncpf  2002).  It  is  nevertheless  difficult  to 
parameterize  the  breaking  of  random  waves,  since  the 
only  available  wave  information  here  are  the  spectral 
densities  l\k,  0).  These  densities  can  be  related  to  the 


orbital  velocity  variance  in  a  narrow  frequency  band. 
This  question  is  addressed  in  detail  by  Filipot  et  al.  (2010) 
Yet.  a  proper  threshold  has  to  be  defined  for  this  quantity, 
and  the  spectral  rate  of  energy  loss  associated  with 
breaking  has  to  be  defined.  Also,  breaking  is  intricately 
related  to  the  complex  nonlinear  evolution  of  the  waves 
(c.g..  Banner  and  Pcirson  2007). 

These  difficulties  will  be  ignored  here.  We  shall  pa¬ 
rameterize  the  spectral  dissipation  rate  directly  from  the 
wave  spectrum,  in  a  way  similar  to  WZB.  Essentially,  we 
distinguish  between  spontaneous  and  induced  breaking, 
the  latter  being  caused  by  large-scale  breakers  over¬ 
taking  shorter  waves,  resulting  in  their  dissipation.  For 
the  spontaneous  breaking  we  parameterize  the  dissipa¬ 
tion  rate  directly  from  the  spectrum,  without  the  inter¬ 
mediate  step  of  estimating  a  breaking  probability. 

We  started  from  the  simplest  possible  dissipation  term 
formulated  in  terms  of  the  direction  integrated  spectral 
saturation  B(k)  given  by  Fq  (5).  with  a  realistic  threshold 
B()r  -  1.2  x  10  1  corresponding  to  the  onset  of  wave 
breaking  (Babanin  and  Young  2005).  This  saturation 
parameter  corresponds  exactly  to  half  of  the  a  paranv 
eter  defined  b\  Phillips  ( 1958).  The  value  B{)  =  4  X  10  \ 
given  by  Phillips,  corresponds  to  a  self-similar  sea  state  in 
which  waves  of  all  scales  have  the  same  shape,  limited 
by  the  breaking  limit. 

This  view  of  the  sea  state,  however,  ignores  com¬ 
pletely  wave  directionality.  Early  tests  of  parameteri 
/aliens  based  on  this  definition  of  B  indicated  that  the 
spectra  were  too  narrow  (Ardhuin  and  Le  Boyer  2006). 
This  effect  could  be  due  to  many  errors.  Because  Banner 
et  al.  (2002)  introduced  a  directional  width  in  their  sal 
uralion  to  explain  some  of  the  variability  in  the  observed 
breaking  probabilities,  we  similarly  modify  the  defini¬ 
tion  of  B.  Expecting  also  to  have  different  dissipation 
rales  in  different  directions,  we  define  a  saturation  B ' 
that  would  correspond,  in  deep  water,  to  a  normalized 
velocity  variance  projected  in  one  direction  (in  the  case 
Sfj  —  2),  with  a  further  restriction  of  the  integration  of 
directions  controlled  by  A„: 


B'(k.e)  = 


-0  »  \ 

k}  cg$s*(6 


C 

o')F(k.ir)-±<w. 

27 T 


(12) 


Here,  we  shall  always  use  A„  SO  .  In  our  model,  a  sea 
state  with  two  systems  with  the  same  energy  but  of  op 
posite  direction  will  thus  produce  less  dissipation  than 
a  sea  state  with  all  the  energy  radiated  in  the  same 
direction. 

Wc  Imally  detine  our  dissipation  term  as  the  sum  of 
the  saturation  based  term  of  Ardhuin  et  al.  (2008a)  and 
a  cumulative  breaking  term  .Shk.ru* 
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r***  ,  max[/J(A)  —  B  ,  01“ 

X  (U)  =  A  '  ' 

[K  [  +  (l  -Srf)max(B'(A\0) 

where 

B(k)-~  max{/T(A\  0).  6  e  [0,  2ir|}.  (14) 

The  combination  of  an  isotropic  part  (the  term  that 
multiplies  6(/)  and  a  direction-dependent  part  (the  term 
with  I  -  8t{)  was  intended  to  allow  some  control  of  the 
directional  spread  in  the  resulting  spectra.  This  aspect  is 
illustrated  in  Fig.  2  with  a  hindeast  of  the  3  November  1999 
ease  during  the  Shoaling  Waves  Experiment  (SHOWEX; 
Ardhuin  et  al.  2(X)7).  Clearly,  the  isotropic  saturation 
in  the  TEST442  dissipation  (with  the  original  threshold 
Br  ~  0.0012)  produces  very  narrow  spectra,  even  though 
it  is  known  that  the  Dl A  parameterization  for  nonlinear 
interactions  tends  to  broaden  the  spectra.  The  same 
pattern  of  behavior  is  obtained  with  the  isotropic  pa¬ 
rameterization  by  van  dcr  Wcsthuysen  et  al.  (2007),  as 
demonstrated  by  Ardhuin  and  Le  Boyer  (2006),  Furl  her, 
using  an  isotropic  dissipation  at  all  frequencies  yields  an 
energy  spectrum  that  decays  faster  toward  high  frequen¬ 
cies  than  the  observed  spectrum  (Fig.  2a).  On  the  con¬ 
trary,  a  fully  directional  dissipation  term  (TEST443  with 
8<t  ~  0)  gives  a  better  lit  for  all  of  the  parameters.  With 
\ti  -  2,  we  reduce  Br  to  0.0009,  a  threshold  for  the  onset 
of  breaking  that  is  consistent  with  the  observations  of 
Banneret  al.  (2000)  and  Banneret  al.  (2002),  as  discussed 
bv  Babanin  and  van  der  Wcsthuysen  (2008).  The  overall 
dissipation  term  in  THST443  is  anisotropic  due  to  the 
cumulative  effect,  but  this  does  not  significantly  alter  the 
underestimation  of  the  directional  spread 

The  dissipation  constant  C^1  was  adjusted  to  2.2  X  ]()  4 
in  order  to  give  acceptable  time-limited  wave  growth 
and  good  directional  fetch  limited  growth  as  described 
by  Ardhuin  et  al.  (2007).  As  noted  in  this  previous  work, 
similar  growth  patterns  in  wave  energy  with  fetch  are 
possible  with  almost  any  magnitude  of  the  wind  input, 
but  a  reasonable  mean  direction  in  slanting  fetch  con¬ 
ditions  selects  the  range  of  possible  lev  els  of  input.  Here, 
the  mean  directions  at  the  observed  peak  frequency  arc 
still  biased  by  about  25°  toward  the  alongshore  direction 
with  the  paramelerizations  proposed  here  (Fig.  2b),  which 
is  still  less  that  the  50c  obtained  with  the  weaker  Tolman 
and  Chalikov  (1996)  source  terms  (Ardhuin  et  al.  2(XI7, 
Fig  1 1 ).  A  relatively  better  ht  is  obtained  with  the  BAJ 
parameterization.  This  is  likely  due  to  either  the  stron¬ 
ger  wind  input  or  the  weaker  dissipation  at  the  peak. 
Both  features  of  the  BAJ  parameterization  are  more 
realistic  than  what  we  propose  here. 


B  ft]:j/7(A'")  +  Vcu(^  *)  +  W*- ">•  03) 


The  equilibrium  sea  state,  achieved  by  most  models 
for  long  durations  with  steady  wind  and  infinite  fetch, 
and  often  compared  with  the  Pierson  and  Moskowit/ 
(1964)  spectrum,  is  largely  controlled  by  the  balance  be¬ 
tween  the  nonlinear  flux  of  energy  to  low  frequencies  and 
the  wind  output  term,  as  discussed  below.  Figure  3  shows 
the  fetch-limited  energy  growth  lor  the  various  param¬ 
elerizations  We  repeal  here  the  sensitivity  lest  to  the 
presence  of  swell  already  displayed  in  Ardhuin  et  al. 
(2(X)7).  Whereas  the  I  m  swell  causes  an  unrealistic 
doubling  of  the  wind  sea  energy  al  short  fetch  in  the  BAJ 
parameterization,  the  new'  parameterizations,  just  like 
the  one  by  van  der  Wcsthuysen  ct  al.  (2007),  are  by 
design  insensitive  to  swell  (not  shown). 

The  dissipation  Sunh  due  to  wave-turbulence  inter¬ 
actions  is  expected  to  be  much  weaker  (Ardhuin  and 
Jenkins  2(X)6)  and  will  he  neglected  here. 

Finally,  following  the  analysis  by  Filipol  et  al.  (2010), 
the  threshold  Br  is  corrected  foi  shallow1  water,  so  that 
BlBr  in  different  water  depths  corresponds  to  the  same 
ratio  of  the  root-mean  square  orbital  velocity  and  phase 
speed.  For  periodic  and  irrotational  waves,  the  orbital 
velocity  increases  much  more  rapidly  than  the  wave  height 
as  it  approaches  the  breaking  limit.  Further,  due  to  non 
linear  distortions  in  the  wave  profile  in  shallow  water,  the 
height  can  be  twice  as  large  as  the  height  of  linear  waves 
with  the  same  energy.  To  express  a  relevant  threshold 
from  the  elevation  variance,  we  consider  the  slope 
k!f\iU(kD)  of  an  hypothetical  linear  wave  that  has  the 
same  energy  as  the  wave  of  maximum  height.  In  deep 
water.  A//lin(x)  ^  ().77,  and  for  other  water  depths  we 
thus  correct  B,  by  a  factor  j  A //|lt>( A*/>)/A //|,n( ^ ) ]  •  Using 
stream  function  theory  (Dalrymple  1974),  a  polynomial 
fit  as  a  function  of  Y  tanh(A/))  gives 

B\  =  BY(  \fAY'  4  M,Y2  4-  M2Y  +  M{ ),  ( IS) 

such  that  B'r  —  Br  in  deep  watci .  The  fitted  constants  are 
\1A  =  1.3286,  V /,  -  -2.3709.  M:  -  1.9995,  and  A/, 
0.2428.  Although  this  behavior  is  consistent  with  the  var¬ 
iation  of  the  depth-limited  breaking  parameter  y  deriv  ed 
empirically  b\  Ruessink  et  al.  (2003).  the  resulting  dissi¬ 
pation  rate  is  not  yet  expected  to  produce  realistic  results 

I  his  value  of  lhe  maximum  equivalent  linear  heighi,  Hhn 
2v  2/-,  wilh  /*.' be  ingi  hoc  leva  lion  variance,  is  smaller  than  i  he  usual 
value  A//  0.88  due  lo  lhe  correction  for  lhe  nonlinear  wave 

profile  for  which  //  v  2/ 
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Pit;.  2.  Wave  spectra  on  3  Nov  1W  al  buoy  X3  (fetch.  39  km: 
wind  speed  U\u  9.4  in  s  ’),  averaged  over  the  time  window 
12(H)- 17(H)  F.S  1 ,  from  observations  and  model  runs,  with  different 
model  piirameteri/ations  (symbols):  BAJ  stands  lor  Bidlot  et  al. 
(2(H).S).  The  (a)  energy,  (b)  mean  direction,  and  (c)  directional  spread 
are  shown.  Phis  figure  is  analogous  to  Figs.  10  and  1 1  in  Ardhuin  et  al 
(2007).  and  the  miKlel  forcing  and  setting  arc  identical.  It  was  further 
verified  that  halving  the  resolution  from  I  km  to  *>00  m  docs  not 
affect  the  results.  All  parametcnzations  settings  arc  listed  in  Tables 
A I  and  A2.  Input  parameters  for  TPS  1443  arc  identical  to  those  for 
I  EST44I,  and  TEST442  differs  from  TES144I  only  in  its  isotropic 
direct  breaking  term,  given  by  %  ft.  180",  and  H,  0.0012. 

for  surf  /ones  because  no  effort  was  made  to  verify  this 
aspect.  This  is  the  topic  of  ongoing  work,  which  is  outside 
of  the  scope  of  the  present  paper. 

The  cumulative  breaking  term  Xhkcu  represents  the 
smoothing  of  the  surface  bv  big  breakers  with  celerity  C 
that  wipe  out  smaller  waves  of  phase  speed  C.  Due  to 
uncertainties  in  the  estimation  of  this  effect  in  the  obser¬ 
vations  of  Young  and  Babanin  (2006),  we  use  the  theo¬ 
retical  model  of  Ardhuin  el  al. (2009b).  Briefly,  the  relative 
velocity  of  the  crests  is  the  norm  of  the  vector  difference, 
A,  ~  C  —  C'|.  and  the  dissipation  rate  of  short  wave  is 
simply  the  rate  of  passage  of  the  large  breaker  over  short 
waves,  that  is,  the  integral  of  A(  A(C)  dC.  where  A(C)  (1C 
is  the  length  of  breaking  crests  per  unit  surface  that  have 
velocity  components  between  Cx  and  Cx  +  dCx%  and 


Non-dimensional  fetch  X* 

FtG.  3.  Fetch  -limited  growth  of  the  vvindsen  energy  ;is  <i  function 
of  fetch  on  3  Nov  1999.  averaged  over  the  time  window  12(H)- 1700 
F.ST,  from  observations  and  model  runs,  with  different  model  pa 
rameterizations  (symbols):  BAJ  stands  for  Bidlot  et  al  (200.S)  This 
figure  is  analogous  to  Fig.  S  in  Ardhuin  et  al  (2007),  the  model 
forcing  and  settings  are  identical.  All  parameters  for  BAJ.TI  S  I  441. 
and  I  ESI  443  arc  listed  in  tables  A I  and  A2  The  input  parameters 
for  TEST443  are  identical  to  those  for  TFST44I,  and  ITS  I  442 
differs  from  I  t  S  1  ‘44 1  only  in  its  isotropic  direct  breaking  term, 
given  bv  v/f  -  0.  lftCT.  and  H  ft.ft()12. 

between  Cy  and  Cv  +  dCx  (Phillips  19<S5).  Because  there 
is  no  consensus  on  the  form  of  \  (Gemtnrich  cl  al.  200S). 
we  prefer  to  link  \  lo  the  breaking  probabilities  Based 
on  Banneret  al.  [(2(HH).  Fig  G ) ; /> 7  -  22(k  -  0.055)2]  and 
taking  their  saturation  panmietci  t  to  be  of  the  order 
of  I.fiy  Br(k.W).  the  breaking  probability  of  dominant 
waves  is  approximately 

P  ~  56.8(max(  \  B'(k.  0)  -  ^  B'r%  0)]2.  ( 16) 

However,  because  they  used  a  zero-crossing  analysis  for 
a  given  wave  scale,  there  are  many  times  when  waves  are 
not  counted  since  the  record  is  dominated  by  another 
scale  in  their  analysis  there  is  only  one  wave  at  any  given 
time.  This  tends  to  overestimate  the  breaking  proba¬ 
bility  by  a  factor  of  1.5-2  (Manasseh  et  al.  2006),  com¬ 
pared  to  the  present  approach  in  which  we  consider  that 
several  waves  (of  different  scales)  may  be  present  at  the 
same  place  and  time.  We  shall  thus  correct  for  this  effeel, 
by  simply  dividing  P  by  2. 

With  this  approach  we  define  the  spectral  density  of 
the  crest  length  (breaking  or  not)  per  unit  surface  /( k) 
such  that  J7(k)  dkx  dkx  is  the  total  length  of  all  crests  per 
unit  surface,  with  a  crest  being  defined  as  a  local  maximum 
of  the  elevation  in  one  horizontal  direction.  In  the  wave¬ 
number  vector  spectral  space  we  take 
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/(k)-  1/(2t7* 2 *A).  (17) 

which  is  equivalent  toa  constant  in  wavenumber-direction 
space  /(A,  0)  =  1  /(2tt2).  This  number  was  obtained  by 
considering  an  ocean  surface  full  of  unidirectional  waves, 
with  one  crest  for  each  wavelength  2  77/A  for  each  spectral 
interval  A k  —  A;  for  example,  one  crest  corresponding  to 
spectral  components  in  the  range  0.5-1. 5  A.  This  potential 
number  of  crests  is  doubled  by  the  directionality  of  the 
sea  state.  These  two  assumptions  have  not  been  verified 
and  thus  the  resulting  value  of  /( k)  is  merely  an  adjustable 
order  of  magnitude. 

Thus,  the  spectral  density  of  the  breaking  crest  length 
per  unit  surface  is  \(k)  =  /(k)P(k).  Assuming  that  any 


breaking  wave  instantly  dissipates  all  the  energy  of  all  of 
the  waves  with  frequencies  higher  bv  a  factor  rCM  or  more, 
then  the  cumulative  dissipation  rate  is  simply  given  by  the 
rate  at  which  these  shorter  waves  are  taken  over  by  larger 
breaking  waves  times  the  spectral  density;  namely. 


hk.cu 


(kJ))  =  CF(k.O) 


r<  r  f 


ArA(k>7k\ 


(18) 


where  /*cu  defines  the  maximum  ratio  of  the  frequencies 
of  long  waves  that  will  wipe  out  the  short  waves. 

We  now  obtain  A  by  extrapolating  Eq.  (16)  to  higher 
frequencies.  With  56.8/(47r)  -  1.44.  it  gives  the  source 
term 


Vcu(*.<>)  =  -1-44C »F(k,6) 


max(  v  B(k\  O’)  -  ./«  . 


We  shall  take  rcu  —  0.5,  and  Ccu  is  a  tuning  coefficient 
expected  to  be  a  negative  number  of  order  1.  w  hich  also 
corrects  for  errors  in  the  estimation  of  /. 

This  generic  form  of  the  source  terms  produces  mark- 
edly  different  balances  for  both  mature  and  fully  devel¬ 
oped  seas. 

Eor  mature  seas,  without  cumulative  effect.  Fig.  4 
shows  that  a  balance  is  possible  that  gives  roughly  the 
same  energy  level  and  w  ind  input  term  up  of  to  0.4  Hz  as 
the  BAJ  parameterization.  However,  the  balance  in  the 
tail  results  in  the  energy  level  decreasing  slower  than/  4 
as  the  dissipation  is  too  weak  compared  to  the  input,  and 
thus  the  nonlinear  energy  (lux  is  reversed,  pumping  civ 
ergy  from  the  tail  to  lower  frequencies. 

The  introduction  of  a  strong  cumulative  term 
(TEST437)  allows  a  balance  at  roughly  the  same  energy 
level.  However,  with  the  present  formulation  this  will 
lead  to  a  dissipation  that  is  too  strong  at  high  frequency 
for  higher  winds.  The  introduction  of  the  sheltering 
effect  via  the  parameter  s„  (details  in  section  2d)  is 
designed  to  reach  a  balance  with  a  weaker  cumulative 
effect. 

The  most  important  qualitative  feature  is  the  lack  of 
a  regular  predefined  shape  for  the  normalized  dissipa¬ 
tion  term.  5oc(A.  0)!E(ki  6).  Whereas  the  shape  given  by 
<3A  +  (1  -  8)k2  is  clearly  visible  in  BAJ  (with  extremely 
high  dissipation  rates  if  one  considers  high  frequencies), 
and  the  low-  to  high-frequency  dissipation  transition  at 
2 fp  is  evident  in  TC.  the  shapes  of  the  new  dissipation 
rates  are  completely  dictated  by  the  local  spectral  satu 
ration  level.  This  leads  to  a  relatively  narrow  peak  of 
dissipation  right  above  the  spectral  peak,  where  the 
saturation  is  strongest. 


This  feature  helps  to  produce  realistic  spectral  shapes 
near  the  peak,  with  a  steeper  low-frequency  side  and 
a  more  gentle  slope  on  the  high-frequency  side,  contrary 
to  the  backward-facing  spectra  produced  by  BAJ  and 
TC.  However,  this  localized  strong  relative  dissipation. 
A'oc(A,  0)/F(k%  0),  is  hard  to  reconcile  with  the  time  and 
spatial  scales  of  breaking  events,  and  is  thus  probably 
exaggerated.  Indeed,  there  should  be  no  significant  dif¬ 
ference  in  the  relative  dissipation  among  the  spectral 
components  that  contribute  to  a  breaking  wave  crest, 
provided  that  they  do  not  disperse  significantly  over  the 
breaker  lifetime,  which  is  less  than  a  wave  period.  There 
is  no  physical  reason  wh>  a  breaking  event  would  take 
much  more  energy,  relatively  speaking,  from  the  spec¬ 
tral  band  ( I .  I  1 .2 )fp  than  from  (1.2  1 .3)/,.  "flic  factor  of 

2  difference  found  here  in  the  relative  dissipation  rates 
thus  appears  to  be  unrealistic.  This  strong  relative  dis¬ 

sipation  at  the  peak  (50%  higher  than  in  BAJ)  is  one 
important  factor  that  leads  to  slower  growth  of  the  wave 
spectrum  in  TEST441  compared  to  BAJ.  It  is  possible 
that  the  localization  of  SIW  at  the  peak  compensates  for 
the  broader  spectrum  produced  by  the  D1A  compared 
to  results  with  an  exact  nonlinear  interaction  calculation. 

We  now  consider  “fully  developed”  conditions,  as  il¬ 
lustrated  by  Fig.  5.  At  low  frequency,  the  nonlinear  swell 
damping  term  Soul  (the  negative  part  of  S,*,,)  cancels 
about  30% -50%  of  the  nonlinear  energy  llux,  so  that  the 
sea  state  grows  only  very  slowly.  As  a  result,  “full  de¬ 
velopment''  does  not  exist,  but  the  resulting  energy  is  still 
compatible  with  the  observations  of  mature  sea  states 
(Alves  et  al.  2003).  In  contrast,  the  linear  swell  damping 
adjusted  by  Tolman  (2(K)2b)  is  designed  to  produce  rea¬ 
sonable  swell  heights  in  the  tropics,  but  it  is  much  smaller 
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spectra  E(f)  and  normalized  source  term  s  M(f)*  S(f) 


0  0.2  0  4  0.6  0.8 


Frequency  (Hz) 


I Ki.  4  Academic  test  case  over  a  uniform  ocean  with  a  uniform 
II)  m  s  wind  starling  from  rest,  after  X  h  of  integration  when 
t'fji I.  Source  term  balances  arc  given  by  the  parameterization 
BAJ,  and  the  parameterizations  proposed  here  with  the  successive 
introduction  of  the  cumulative  breaking  and  the  wind-sheltering 
effects  with  the  parameters  Ccn  and  s„.  For  BA.I.  a  diagnostic/  * 
tail  is  applied  above  2.5  x  the  mean  frequency.  The  source  terms 
are  multiplied  by  lhe  normalization  function  \1(f)  p„C7 
[()„£(  f)trUu\].  The  TC  stands  for  of  Tolman  and  Clmlikov  ( 1 996 ) . 

than  the  nonlinear  energy  flux  to  low  frequencies,  even 
with  the  reduced  interaction  coefficient  proposed  by 
Tolman  and  Chalikov  (1996).  A  nonlinear  swell  dissipa¬ 
tion  appears  to  be  necessary  to  obtain  both  a  realistic 
damping  of  observed  swells  and  a  satisfactory  agreement 


"dr  cto  da  dr 


Ccu—0,  su“^ 

M . dv“d8 


T15 


Ccu  su=0 

~di — d2 — da — to — d5 — 07 — d8 


Q5 

0 

-Q5 


Ccu  —0-4,  Su=l 


FlG.  5.  As  in  Fig.  4.  but  after  4K  h  of  integration  and  without 
normalization  of  the  source  terms.  Source  term  balances  are  given 
bv  the  parameterization  BAJ  and  the  parameleri'/ations  proposed 
here  with  the  successive  introduction  of  the  cumulative  breaking 
and  the  wind  sheltering  effects  with  lhe  parameters  (  and  s„ 


with  mature  wind  waves.  Nonlinearity  also  brings  within 
the  same  order  of  magnitude  the  decay  scales  estimated 
for  short  (Hogstiom  et  al  2009)  and  very  long  swells 
(Ardhuin  et  al.  2009b). 

Both  parameten/ations  are  physically  very  different 
from  the  parameterization*  of  the  KHH  family,  includ¬ 
ing  Bidlot  et  al.  (2005).  In  these  eases,  the  swell  energy  is 
tost  to  the  ocean  via  whitecapping.  Here,  we  propose 
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that  this  energy  is  lost  to  the  atmosphere,  with  an  asso¬ 
ciated  momentum  flux  that  drives  the  wave-driven  wind 
observed  in  wave  flumes  (Harris  1966)  and  for  very  weak 
winds  at  sea  (Smedman  et  al.  2009). 

In  the  inertial  range,  a  reasonable  balance  of  all  the 
source  terms  is  obtained  for  Ccu  =  -0.4  (Fig.  5).  In  this 
ease,  the  spectrum  approaches  an /  4  shape  up  to  0  4  Hz. 
The  behavior  of  the  high-frequency  tail  is  best  seen  when 
displayed  in  nondimensional  form,  as  is  done  in  Fig.  6. 
Figure  6  shows  the  unrealistically  high  level  of  the  tail 
without  cumulative  effect  nor  modification  of  the  wind 
input,  for  fully  developed  waves.  With  the  TEST441  pa¬ 
rameterization  (Ccu  ~  0.4  and  .s„  =  I),  the  spectral  level  al 
0.7  Hz  (3  m  wavelength)  may  be  30%  loo  high  com¬ 
pared  to  other  analyses  by  Long  and  Resio  (2007),  but 
the>  are  in  the  range  0.0045-0.008  given  by  Babanin  and 
Soloviev  (1998).  Further,  for  average  oceanic  condi¬ 
tions,  the  resulting  spectral  moments  m 3  is  still  slightly 
underpredicled  (Ardhuin  ct  al.  2009b).  Spectral  levels 
for  shorter  waves  (20  cm-2  m)  were  estimated  at  a  ~ 
1 .8  X  10  3  by  Banner  ct  al.  (1989).  There  is  thus  a  strong 
need  for  more  spectral  measurements  in  the  range  of 
wavelengths  from  0.5  to  5  m,  outside  of  the  range  of  buoy 
measurements. 

Adding  the  cumulative  effect  can  be  a  means  of  con¬ 
trolling  the  tail  level,  but  this  degree  of  freedom  is  not 
enough.  Indeed,  in  strongly  forced  conditions  the  dom¬ 
inant  waves  break  frequently,  and  a  high  cumulative 
effect,  Ccll  =  —  1 ,  reduces  the  energy  level  in  the  tail  below 
observed  levels.  This  effect  can  be  verified  with  mean 
square  slopes  estimated  from  satellite  altimeter  mea¬ 
surements  (Fig.  8),  or  high  moments  of  the  frequency 
spectrum  derived  from  buoy  data  (not  shown  hut  similar). 

A  more  accurate  shape  of  the  spectrum  tail  may  be 
obtained  with  a  lower  Ccl)  or  a  higher  rtu.  so  that  domi¬ 
nant  breaking  waves  will  only  wipe  out  much  smaller 
waves.  Instead,  and  because  the  wind  to  wave  momentum 
flux  was  loo  high  in  high  w  inds,  we  chose  to  introduce  one 
more  degree  of  freedom,  allow  ing  a  reduction  of  the  w  ind 
input  at  high  frequency. 

d.  Wind  in  pm 

The  wind  input  parameterization  is  thus  adapted  from 
Janssen  ( 1991;  sec  also  Chalikov  1993)  and  the  following 
adjustments  performed  by  Bidlot  et  al  (2005.  2007a). 
The  full  wind  input  source  term  reads 


9)  - 


x  max  [cos  (ft  —  0M),  oyW-’fA,  ft). 


(19) 


where  /3m;lv  is  a  nondimensional  growth  parameter  (con¬ 
stant).  tt+  is  the  wind  friction  velocity,  and  k  is  von 


Fus.  6.  Values  of  the  speciral  sutural  ion  for  the  cases  presented 
in  Fig.  5. 


Kami  an  *s  constant.  In  the  present  implementation  the 
air-water  density  ratio  is  constant  The  power  of  the 
cosine  is  taken  constant  with  p  2.  We  define  the  ef¬ 
fective  wave  age  /  logt/x),  vv  here  p  is  given  by  Janssen 
and  corrected  for  intermediate  water  depths,  so  that 

Z  -  log(Az,)  +  K/[cos(ft  -  0u)(uJC  +  zjl  (20) 

where  Z\  is  a  roughness  length  modi  lied  by  the  wave- 
supported  stress  r„  and  :  ,  is  a  wave  age  tuning  param¬ 
eter.  We  define  z\  implicitly  b\ 

7 '“<:*)•  l2l) 


*11  '  '■()  m.iv 


and 


(22) 


*0 


(23) 


where  is  the  height  at  which  the  wind  speed  is  speci¬ 
fied,  usually  10  m.  The  maximum  value  of  Co  was  added 
to  reduce  possible  unrealistic  w  ind  stresses  al  high  winds 
that  are  otherwise  given  by  the  standard  parameteriza¬ 
tion.  For  example,  cUmilx  —  0.0015  is  equivalent  to  set¬ 
ting  a  maximum  wind  drag  coefficient  of  2.5  X  10  \  For 
the  TEST44I  parameterization,  we  have  adjusted  zn 
0.006  and  =  1.52  (Fig.  7). 

An  important  part  of  the  parameterization  is  the  cal¬ 
culation  of  the  w  ave  supported  stress  rH.,  which  includes 
the  resolved  part  of  the  spectrum,  as  well  as  the  growth 
of  an  assumed  /  ^  diagnostic  tail  beyond  the  highest 
frequency.  This  parameterization  is  highly  sensitive  to  the 
high-frequency  part  of  the  specl mm  since  a  high  energy 
level  there  will  lead  to  a  larger  value  of  z i  and  //*,  which 
gives  a  positive  feedback  and  reinforces  the  energy  levels. 
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R(i.  7.  lncrcmcnt.il  adjustments  lo  the  wind  wave  interaction 
source  term  .Salm,  going  from  the  13  AJ  form  to  the  one  used  in 
I I4S 1 44 1 .  Errom  one  curve  lo  the  ncxl  only  one  parameter  is 
changed.  Differeni  are  given  for  the  same  wave  spectrum, 
previously  obtained  by  running  the  model  from  a  calm  sea  for  K  h 
with  ihe  13 A  )  parameterization  and  a  wind  speed  of  10  m  s  r. 


To  allow  a  balance  with  the  saturation  based  dissi 
pation,  the  wind  input  at  high  frequency  is  reduced  by 
modifying  the  friction  velocity  //*.  This  correction  also 
reduces  the  drag  coefficient  at  high  winds.  In  particular, 
the  wind  input  is  reduced  for  high  frequencies  and  high 
winds,  loosely  following  Chen  and  Belcher  (2000).  This 
is  performed  by  replacing  t/*  in  Eq.  (19)  with  a  wave- 
number  dependent  */*(A‘)  defined  bv 


(//*):  =  ni(coR0|#,  sin#M) 

f*r2’-u*'.0> 

r  <r  I 


(cos 0.  fine) ilk1  <10.  (24) 


where  the  sheltering  coefficient  .v„|  —  1  can  be  used  lo 
tunc  the  wind  stresses  for  high  winds,  which  would  be 
largely  overestimated  for  sn  —  0.  For  sn  >  0,  this  shel¬ 
tering  is  also  applied  within  the  diagnostic  tail,  which 
requires  the  estimation  of  a  three-dimensional  lookup 
table  for  the  high-frequency  stress.  The  shape  of  the  new 
w  ind  input  is  illustrated  in  Fig.  7  for  fully  developed  seas. 
Clearly,  for  relatively  young  waves  the  energy  levels  at 
the  spectral  peak  are  lower  with  sf<  =  1  (TEST441 )  than 
in  other  runs;  this  is  largely  clue  to  a  reduced  feedback  of 
the  wave  age  on  the  wind  stress  via  the  tJt  term  in  Fq. 
(23).  fhe  reduction  of  **:<,  from  0.01 1  lo  0.006  strongly 
reduces  the  input  for  frequencies  in  the  range  0.15- 
0.2  Hz,  which  is  probably  overestimated  in  BAJ  when 
average  (20%)  levels  of  guslincss  are  considered:  now 
l ho  wind  input  goes  lo  zero  for  /  =  0.13,  which  corre¬ 
sponds  to  CHI  in  =  0.83,  whereas  it  is  still  significant  at 


that  wave  age  in  the  BAJ  parameterization.  As  a  result 
the  much  lower  input  level  needs  a  readjustment,  per¬ 
formed  here  by  increasing  pmux  lo  1.52.  Yet,  this  high 
value  of  /3mav  produces  very  high  wind  stress  values  and 
thus  a  very  strong  high-frequency  input.  Adding  the 
sheltering  term  s„  ~  1  allows  a  decent  balance  at  high 
frequency.  Finally,  the  addition  of  the  air-sea  friction 
term  that  gives  swell  dissipation  produces  a  significant 
reduction  of  the  input  lo  the  wind  sea  at  f  =  0.25  Hz.  It  is 
questionable  w  hether  this  mechanism  also  applies  in  the 
presence  of  the  critical  lay  er  for  those  waves.  This  mailer 
clearly  requires  more  theoretical  and  experimental 
investigation. 


3.  Consequences  of  the  source  term  shape 

We  have  already  illustrated  the  effects  of  various  pa 
ramelcrs  on  spectral  shapes  in  academic  lime-limited 
and  more  realistic  fetch-limited  conditions.  We  now  look 
at  real  sea  stales  observed  in  the  World  Ocean.  Although 
wave  spectra  arc  difficult  to  compare  to  the  few  available 
observations,  we  have  investigated  the  systematic  varia¬ 
tion  of  spectral  moments 


r /•  (/.  9)d9df, 


(25) 


with  n  -  2,  3,  and  4,  and  cutoff  frequencies  in  the  range 
0.2-0.4  Hz.  The  spectral  density  F(j\  0)  is  estimated  as 
2 it {{ k .  0)iC^  using  linear  wave  theory.  Such  moments 
arc  relevant  to  a  variety  of  applications.  Ardhuin  cl  al. 
(2(X)9b)  investigated  the  third  moment,  which  is  pro 
porlional  lo  the  surface  Stokes  drill  in  deep  water,  and 
found  that  buoy  data  are  very  well  represented  by  a  sim¬ 
ple  function  of  U m.  //  and  which  typically  explains 
95%  of  the  variance  of  nu. 

TTiis  relationship  is  well  reproduced  in  hindcasls  using 
Ccu  0.4  and  .v,f  =  f ,  while  the  BAJ  source  terms  give 
almost  a  constant  value  of  //n  when  //,  varies  and  U in  is 
fixed  (Ardhuin  el  al.  2009b). 

Here,  we  consider  the  fourth  moment  w4.  which,  for 
linear  waves,  is  proportional  lo  a  surface  mean  square 
slope  filtered  al  the  frequency  j\  .  For  modeled  values 
a  constant  0.011  is  added  lo  account  for  the  short  waves 
that  contribute  lo  the  satellite  signal  and  that  arc  not  re¬ 
solved  in  the  model.  This  saturated  high  frequency  tail, 
independent  of  wind  speed,  is  consistent  with  the  ob¬ 
servations  of  Banneret  al.  ( 1989)  and  Vandemark  el  af 
(2004). 

Figure  8  shows  that  lor  any  given  wind  speed  mss< 
increases  with  the  wa\c  height  (Gourrion  el  al.  2002). 
whereas  this  is  not  the  case  for  m4  in  the  BAJ  parame¬ 
terization.  or.  for  ver\  high  winds,  when  Ccu  is  too  strong. 
In  the  case  of  BAJ,  this  is  due  lo  the  (klk,Y  part  in  the 
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El( i.  8.  Variation  of  the  surface  mean  square  slope  estimated  as  either  (bottom  left)  0,64  <rn  using  the  C -hand 
altimeter  on  board  Jason  I  after  the  correction  of  a  1.2-dB  bias  in  the  Jason  data,  or  bv  integration  of  modeled  spectra 
from  0  to  0.72  117.  with  either  (top  left)  the  new  parameterization  with  Ctu  1  and  \  I.  (bottom  right)  with 
C\u  0.4  and  st(  I  (TLST44I )  or  (top  right)  the  parameterization  by  BAJ.  The  original  1  \\7  data  from  Jason  is 
subsampled  at  0.5  H/  and  averaged  over  10  s  (i.e.,  58  km  along  the  satellite-  track).!  he  same  averaging  is  applied  to 
the  wave  model  result,  giving  the  total  of  393  3S2  observations  reported  here,  for  the  first  half  year  of  2007. 


dissipation  term  [Eq.  (2)],  which  plays  a  rale  similar  to 
the  cumulative  term  in  our  formulation.  For  Ccll  =  -1 
and  v,,  —  0.  the  cumulative  effect  gets  too  strong  for  wind 
speeds  over  10  m  s  \  in  which  case  starts  to  decrease 
with  increasing  wave  height,  whereas  for  high  winds  and 
low  (t.e.,  young)  waves,  the  high-frequency  tail  is  too 
high  and  the  mean  square  slope  gets  to  be  as  large  as  6%. 
which  is  unrealistic.  For  -  0  the  high-frequency  tail 
responds  too  much  to  the  wind,  hence  our  use  of.vM  =  1  in 
the  TEST44I  combination.  The  presence  of  a  cuniulativ e 


dissipation  term  allows  for  a  different  balance  in  the 
spectral  regions  above  the  peak,  where  an  equilibrium 
range  with  a  spectrum  proportional  to/  4  develops  (l  ong 
and  Resto  2007),  and  in  the  high-frequency  tail  were  the 
spectrum  decays  like/  or  possibly  a  little  faster.  The 
spectral  level  in  the  range  0.2-0.4  Hz  was  carefully  com¬ 
pared  against  buoy  data  and  was  found  to  be  realistic. 

These  interpretations  relv  on  the  assumption  that  the 
model  results  are  really  wavenumber  spectrum  E(k%  0) 
converted  to  frequency  densities  using  linear  wave 
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theory,  which  is  expected  to  be  closer  to  the  frequency 
spectra  obtained  from  Lagrangian  buoy  measurements, 
and  the  wavenumber  spectrum  of  nonlinear  waves 
(Janssen  2(X)9).  This  matter  is  left  for  further  studies,  to¬ 
gether  with  a  detailed  interpretation  of  the  altimeter  radar 
cross  sections.  Although  it  covers  much  less  data,  the 
analysis  of  m4  obtained  from  buoy  heave  spectra  pro¬ 
duces  results  similar  to  Fig.  8. 


4.  Verification 

To  provide  simplified  measures  of  the  difference  be¬ 
tween  the  model  time  scries  A'ni<Ki  and  the  observations 
A'obs,  we  use  the  following  definitions  for  the  normalized 
root- mean  square  error  (NRMSK), 


NRMSH(A')  = 


(26) 


the  normalized  bias. 


NB(A') 


ohs  ^  nunl  , 


and 


(27) 


Pearson's  linear  correlation  coefficient. 


r(X)  - 


ohs  obs  ^  ^  m<Hl 


mix!  mod  ) 


(28) 


where  the  overbar  denotes  the  arithmetic  average. 

The  normalization  of  the  errors  allows  a  quantitative 
comparison  between  widely  different  sea  state  regimes. 
Because  previous  studies  have  often  used  (nonnornializcd) 
RMSEs.  we  also  provide  RMSE  values.  In  addition  to  the 
coastal  fetch  limited  case  of  SHOWKX,  presented  above, 
the  paramcleri/alions  are  calibrated  on  at  the  global  scale 
and  validated  in  two  other  cases. 

a.  Global-scale  results 

We  present  here  results  for  the  entire  year  2007.  using 
a  stand-alone  0.5°  resolution  grid,  covering  the  globe 
from  80°  south  to  80°  north.  The  model  has  actually  been 
adjusted  to  perform  well  over  this  dataset,  but  the  very 
large  number  of  observations  (more  than  2  million  al¬ 
timeter  collocation  points)  makes  the  model  robust  and 
an  independent  validation  on  2008-10  gives  identical 
results.  The  interested  reader  may  also  look  at  the 
monthly  reports  for  the  Service  Hydrographique  et 
Oceanographique  de  la  Marine  (SHOM)  model  (e.g.. 


Bidlot  2(>08) ,  generated  as  part  of  the  model  verifica¬ 
tion  project  of  the  Intergovernmental  Oceanographic 
Commission-World  Meteorological  Organization  (!OC- 
WMO)  Joint  Commission  on  Oceanography  and  marine 
Meteorology  (JCOMM).  in  which  the  TEST44I  pa¬ 
rameterization  (Ccu  ~  0.4  and  .v„  -  1)  is  used,  except 

for  the  Mediterranean  where,  the  TEST405  has  been 
preferred  for  its  superior  performance  for  younger  seas. 
These  SHOM  models  are  run  in  a  combination  of  two- 
way  nested  grids  (Tolman  2007).  1  he  monthly  JCOMM 
reports  include  both  analyses  and  forecasts,  and  many 
SHOM  forecasts  from  December  2(X)S  to  June  2(X)9  have 
been  affected  by  wind  file  transfer  problems. 

Comparing  the  model  results  for  //v  to  well  calibrated 
altimeter  derived  measurements  (Qucffeulou  and  Croi/.c- 
Fillon  2008)  provides  a  good  method  of  verification  for 
the  model  performance  in  a  number  of  different  wave 
climates.  Figure  9  shows  that,  as  expected,  the  impor¬ 
tant  positive  bias  in  the  swell-dominated  regions  when 
using  the  BAJ  parameterization  has  been  largely  re¬ 
moved.  This  is  essentially  the  signature  of  the  specific 
swell  dissipation  that  is  parameterized  in  The 

largest  bias  pattern  now  appears  in  the  Southern  Ocean, 
reaching  30  cm  in  the  southern  Atlantic  in  2(X)7  and 
more  than  80  cm  in  2004.  Although  this  bias  is  small 
compared  to  the  local  averaged  wave  height,  it  is  rather 
odd  when  the  model  errors  are  plotted  as  a  function  of 
wave  height  in  Fig.  10.  Win  would  the  model  over¬ 
estimate  the  Southern  Ocean  waves  but  underestimate 
the  very  large  waves? 

The  structure  of  the  large  bias,  also  seen  in  model  re¬ 
sults  of  BAJ.  is  reminiscent  of  the  observed  pattern  in 
iceberg  distribution  produced  by  T burnadre  et  al  ( 2( K >8 ) 
for  the  vear  2004.  These  observed  iceberg  distributions 
are  enough  to  give  a  cross  section  for  incoming  waves  of 
the  order  of  1%  10%  for  a  250-km  propagation  length. 
Assuming  a  full  blocking  of  wave  energy  by  icebergs 
larger  than  I  km  in  diameter  is  enough  to  cancel  most  of 
the  bias. 

Also  noticeable  is  a  significant  negative  bias  in  the 
equatorial  South  Pacific,  amplified  from  the  same  bias 
obtained  with  the  BAJ  parameterization.  It  is  possible 
that  the  masking  of  subgnd  islands  (Tolman  2(X)3)  in¬ 
troduces  a  bias  by  neglecting  shoreline  reflections.  This 
model  defect  could  be  exacerbated  in  this  region  by  the 
very  large  ratio  of  shoreline  length  to  sea  area.  This  will 
also  require  further  investigation.  Finally,  the  negative 
biases  for  Hs  along  midlatitude  cast  coasts  arc  reduced 
hut  still  persist.  It  is  well  known  that  these  areas  arc  also 
characterized  by  strong  boundary  currents  (such  as  the 
Gulf  Stream,  Kuroshio.  and  Agulhas  current)  with  warm 
waters  that  arc  generally  conducive  to  wave  amplifica¬ 
tion  and  faster  wind-wave  growth  (e.g.,  Vandemark  el  al. 
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I  k;.  9.  The  bias  (cm)  for  ihc  year  2007.  The  global  0.5  WWATCH  model  is  compared  In 
allimctcrs  aboard  the  Jason -/,  the  Environmental  Satellite  (Envisat).  and  Gcosai  f  ollow  On 
(GFO)  sal  ell  ii  c  following  Ihe  mclhod  of  Rasclc  cl  al.  (2008).  (top)  The  rcsuli  with  ihc  BAJ 
parameterization,  and  (holiom)  Ihc  result  with  the  Ccu  0.4  and  s„  I  (THST441) 
parameterization 


2001 ).  Neither  effect  is  included  in  the  present  calculation 
because  the  accuracy  of  both  modeled  surface  currents 
and  air-sea  stability  parameterizations  arc  likely  to  be 
insufficient  (Collard  et  al.  2008;  Ardhuin  et  al.  2007). 

I  he  reduction  of  systematic  biases  clearly  contributes 
to  the  reduction  of  RMSEs,  as  is  evident  in  the  equa¬ 
torial  eastern  Pacific  (Fig.  11).  However,  the  new  pa¬ 
rameterization  also  brings  a  considerable  reduction  to 
the  scatter,  with  reduced  errors  even  where  biases  are 
minimal,  such  as  the  trade  winds  area  south  of  Hawaii, 
where  the  NRMSE  for  Ifs  can  be  as  low  as  5%.  When 
areas  within  400  km  from  continents  are  excluded,  be 
cause  the  global  model  resolution  may  be  inadequate, 
significant  errors  (>12.5%.  in  yellow  to  red)  remain  in 
the  northern  Indian  Ocean,  on  the  North  American  and 
Asian  cast  coasts,  and  in  the  southern  Atlantic.  The 
parameterizations  TEST405,  TEST437,  and  TEST44I 
produce  smaller  errors  on  average  than  BAJ.  It  is  likely 
that  the  model  benefits  from  the  absence  of  swell  in¬ 
fluences  on  wind  seas:  swell  in  BAJ  typically  leads  to 
a  reduced  dissipation  and  stronger  wind  wave  growth. 
As  models  are  adjusted  to  average  sea  state  conditions, 
this  adjustment  leads  to  reduced  wind  sea  growth  along 
east  coasts,  where  there  is  generally  less  swell. 

For  the  highest  waves,  the  model  scatter  is  smallest, 
with  some  important  biases  for  some  parameterizations. 


Although  the  altimeter  estimates  are  not  expected  to  be 
valid  for  //v  larger  than  about  12  m.  due  to  the  low 
backscatter  and  a  waveform  that  is  too  short  to  properly 
estimate  //,,  this  good  model  performance  contradicts 
some  claims  (e.g.,  CavaJcri  2009)  that  models  may  he 
less  accurate  for  severe  weather. 

Although  much  more  sparse  than  the  altimeter  data, 
the  in  situ  measurements  collected  and  exchanged  as 
part  as  the  JCOMM  wave  model  verification  are  very 
useful  for  constraining  other  aspects  of  the  sea  state 
This  is  illustrated  here  with  mean  periods  Tiu{)z  for  data 
provided  by  the  U.K.  and  French  meteorological  ser¬ 
vices,  and  peak  periods  I'f,  for  all  other  sources.  It  is 
worth  noting  that  the  errors  in  //v  for  in  situ  platforms 
are  comparable  to  the  errors  found  for  altimeter  data. 

With  the  BAJ  parameterizations.  the  largest  errors  in 
the  model  results  are  the  1 .2-1 .8  s  biases  in  peak  periods 
along  the  IJ.S.  west  coast  (Fig.  12),  and  the  under¬ 
estimation  of  peak  periods  along  the  U.S.  east  coast. 
However,  peak  and  mean  periods  off  the  European 
coasts  are  generally  very  well  predicted.  When  using  the 
TEST44I  parameterization,  the  explicit  swell  dissipa¬ 
tion  reduces  the  bias  over  periods  along  the  U.S.  west 
coast,  but  the  problem  is  not  completely  solved,  with 
residual  biases  of  0.2-0.4  s.  This  is  consistent  with  the 
validation  using  satellite  s\  nthetic  aperture  radar  (SAR) 
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Eui  10  Wave  model  errors  as  a  funciion  of  observed  //,:  (top  to 
bottom)  bias,  scatter  index  (St),  NRMSE,  and  number  of  obser¬ 
vations.  The  model  output  ill  3  h  intervals  is  compared  to  that  from 
Jason- /.  I  nvisat.  and  Of  O  following  the  method  of  Rascle  cm  al 
(2008).  Namely,  the  altimeter  I -Hz  Ku  band  estimates  of  //,  arc 
averaged  over  l  This  provides  2  044  54S  observations 


data  (Fig.  I ).  which  showed  a  tendency  to  underpredict 
steep  swells  near  the  storms  and  overprediet  them  in  the 
far  field.  A  simple  increase  in  ihe  swell  dissipation  was 
tested  but  it  tended  to  deteriorate  the  results  of  other 
parameters;  this  trend  is  possibly  associated  with  an 
underestimation  of  wave  heights  in  ihe  center  of  severe 
storms  using  TEST441 .  Along  European  coasts,  despite 
a  stronger  bias.  Ihe  errors  in  Tm{l2  are  particularly  re¬ 
duced,  Again,  this  reduction  of  the  model  scatter  can  be 
largely  attributed  to  ihe  decoupling  of  swell  from 
wmdsea  growth. 

The  general  performance  of  the  parameterizalions  is 
synthesized  in  Table  1.  It  is  interesting  to  note  that  ihe 
parameterization  TEST405,  which  uses  a  diagnostic  l ail 
for  2.3  limes  the  mean  frequency,  gives  good  results  in 
terms  of  the  scatter  and  bias  even  for  parameters  related 
to  short  waves  {niy%  mA).  I  his  use  of  a  diagnostic  tail  is 
thus  a  reasonable  pragmatic  alternative  to  the  more 
costly  explicit  resolution  of  shorter  waves,  which  re¬ 
quires  a  smaller  adaptive  time  step,  and  more  complex 
parameterizalions.  The  diagnostic  tail  generally  mimics 
the  effects  of  both  the  cumulative  and  sheltering  effects. 
Yet,  the  parameterization  TEST44I  demonstrates  that 
it  is  possible  to  obtain  slightly  better  results  with  a  free 
tail.  The  normalized  biases  indicated  for  the  mean 
square  slopes  are  only  relative  because  of  the  approxi¬ 
mate  calibration  of  the  radar  cross  section.  They  show 
that  the  BAJ  parameterization  (Bidlot  ct  al.  2003),  and 
to  a  lesser  extent  the  use  of  a  f  s  tail,  produce  energy 
levels  that  are  relatively  lower  al  high  frequency. 

b.  Lake  Michigan 

At  the  global  scale,  the  sea  stale  is  never  very  young, 
and  il  is  desirable  to  also  verify  the  robustness  of  the 
parameterization  in  conditions  that  arc  more  represen 
lalixe  of  the  coastal  ocean  We  thus  follow  ihe  analysis 
of  ihe  performance  of  wave  models  in  Rogers  and  Wang 
(2007),  hereinafter  RW2007,  and  give  results  for  Lake 
Michigan,  which  is  representative  of  relatively  young 
waves.  The  model  was  applied  with  the  BAJ,  TC. 
TEST437,  and  TEST441  paramcterizations  over  the  same 
time  frame  as  was  investigated  by  RW2007: 1  September 
14  November  2002.  The  model  is  run  with  4-km  resolution 
and  with  other  sellings  and  forcing  fields  similar  to  ihose 
defined  by  RW2007  (c.g..  10°  directional  resolution  and 
wind  field  derived  from  in  siiu  observations) 

Using  the  directional  validation  method  proposed  by 
these  authors,  the  TC  parameterization  underestimates 
the  directional  spread  crH  by  1.2°-1.6°  in  the  range  (0.8 
2.0)//M  and  more  at  higher  frequencies.  The  underes¬ 
timation  with  BAJ  is  about  half,  and  the  TEST441  and 
TEST437  simulations  ovcrprcdicl  the  directional  spread 
by  about  2.3°-5.9°  in  the  range  (f)  8-2.0)/,,,  and  less  so  for 
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Fin.  11.  Normalized  RMSE  (%)  for  the  significant  wave  height  over  the  year  2007  The 
global  (1 5  resolution  WAV  A  I’CH  model  is  compared  to  altimeters  aboard  the  Jason  /.  fan  isat , 
and  (jbO  following  the  method  of  Rascle  el  al.  (2008).  (lop)  Fhe  result  with  lilt*  RA.I  parame 
terization,  and  (bottom)  the  result  with  the  (\lt  =  0.4  and  s„  T  1  ( I  BS  1441 )  parameterization 


higher  Frequencies.  It  thus  appears  that  the  broadening 
introduced  to  fit  the  SHOWEX  1.999  observations  may 
not  be  optimal  for  other  situations. 

Further  results  are  presented  in  Table  2  and  Fig.  13. 
The  top  panel  in  Fig.  1 1  shows  a  comparison  between 
the  summed  values  of  collocated  model  and  observed 
spectral  densities  for  the  duration  of  the  simulation. 
This  presentation  provides  a  frequency  distribution  of 
the  bias  of  the  various  models,  while  also  indicating  the 
relative  contribution  of  each  frequency  to  the  wave  cli¬ 
mate  for  this  region  and  time  period.  The  bottom  panel 
in  Fig.  13  shows  the  correlation  coefficient  r  for  the 
equivalent  significant  wave  heights  computed  for  mill 
tiple  frequency  bands.  This  is  presented  in  terms  of ///;, 
(bin  width  =  0.1 ).  with  fp  being  calculated  as  the  stabi¬ 
lized  “synthetic  peak  frequency"  of  the  corresponding 
buoy  spectrum,  as  defined  in  RW2U07. 

The  most  noticeable  outcome  of  these  comparisons  is  the 
relatively  poor  performance  of  the  TC  paramelcri/ations 


I  aken  in  context  with  other  TC*  results  presented  herein 
and  a  prior  undocumented  application  of  the  model  in 
the  Great  Lakes  with  model  wind  fields,  this  suggests  that 
these  paramelcri/ations  have  sonic  undesirable  depen¬ 
dence  on  scale,  with  the  parameters  adjusted  by  Tolinan 
(2002b)  being  most  optimal  for  ocean  scale  applications 
and  nonoptimal  for  small  scale  applications. 

The  KHH  dissipation  scheme  performs  well.  This  is 
consistent  with  prior  published  applications  with  the 
Simulating  Waves  Nearshore  (SWAN)  model,  (Rogers 
et  al.  2003)  and  RW2007  However,  as  noted  by  those 
authors,  the  KHH  physics  have  an  advantage  in  this 
simple  wave  climate,  and  should  not  be  expected  to 
perform  well  in  mixed  sea  swell  conditions 

The  BAJ,  TEST437  and  I  KST441  models  also  perform 
well  here.  Taken  together  with  the  global  comparisons 
above,  wc  observe  no  apparent  dependency  of  model  skill 
on  the  scale  of  the  application  with  these  three  physics. 
In  the  bias  comparison  (Fig.  13,  top),  the  BAJ  model  is 
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Fl(i  12  Statistics  for  the  year  2007  based  on  the  JCOMM  verification  database  (Bid lot  et  al 
2007b)  Shown  are  the  (a).(d)  bias  and  (h),(c)  NRMST  tor  / p  or  at  in  situ  locations  using 
(a).(h)  the  RAJ  method  or  (c).(d)  the  proposed  TTS  I  441  (  / r  is  shown  at  all  buoys  except  tor 
the  l  J.K  and  French  buoys,  for  which  /mn:  is  shown)  Fhc  different  symbols  are  only  used  to 
help  distinguish  the  various  colors  and  do  not  carry  extra  information. 


nearly  identical  to  the  KHH  model  Similarly,  the  two 
new  models  are  nearly  identical.  In  terms  of  lts  (  Table  2), 
the  TEST441  and  TEST437  yield  minor  underest  ima- 
tions,  but  give  slightly  more  accurate  results  compared 
to  observations  than  BAJ  or  KHH. 

It  thus  appears  that  for  such  young  seas,  the  directional 
spreading  of  the  parameterization  could  be  improved,  but 
the  energy  content  of  various  frequency  bands,  and  as 


a  result  ihc  mean  period,  arc  reproduced  with  less  scatter 
than  with  previous  parameleri/ations 

c.  Hurricane  Ivan 

Although  the  global  hindcasl  does  contain  quite  a  few 
extreme  events,  with  significant  wave  heights  up  to  17  m, 
these  were  obtained  with  a  relatively  coarse  wave  model 
grid  and  wind  forcing  (0.5  resolution  and  6-h  time  step) 
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I  Alii. I*.  I.  The  model  accuracy  for  significant  wave  height  mea¬ 
sured  wave  parameters  over  the  oceans  in  2007.  The  mss  data  from 
Jason- 1  corresponds  to  January  July  2007  (393  3S2  collocated 
points),  while  H ,  statistics  are  obtained  for  the  entire  year  and  all 
available  satellites.  These  global  averages  are  area  weighted,  and 
the  scatter  index  (SI)  and  NRMSE  are  the  area-weighted  averages 
of  the  local  Si  and  NRMSE.  The  Tr  and  Tmo?  statistics  are  averages 
of  statistics  computed  for  each  buoy  separately.  For  we  used 
WMO  buoys  41002. 4 10 10. 42001. 42002, 42003. 44004.  44008.  44011. 
44 1 37.  44 1 38. 44 1 39.  44141.46001. 46004.  46035.  46066.  46 1 84. 46002. 
46005.  46036.  46059.  51001.  51002.  51003.  and  51004;  for  7m(C.  we 
used  62029.  62081. 62163.  and  64045;  and  for  niu  we  Ux>k  the  results 
tor  buoy  46(X)5  as  given  in  Ardhuin  el  al. (2009a).  Unless  otherwise 
specified  by  the  number  in  parenthesis  the  cutoff  frequency  is  taken 
to  be  0.4  Hz,  where  C  is  for  C  band  The  normalized  bias  (NO)  is 
defined  as  the  bias  di\  ided  by  the  RMS  observed  value,  while  the  SI  is 
defined  as  the  RMS  difference  between  the  modeled  and  observed 
values,  after  correct  ion  for  the  bias,  normalized  by  the  RMS  observed 
value,  and  r  is  Pearson's  correlation  coefficient 
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T  \m  l-  2.  Model  data  comparison  ai  NDBC  buoy  45007  (Lake 
Michigan)  for  four  frequency  bands  Statistics  are  given  for  equiv¬ 
alent  significant  wave  heights  in  the  bands  0.5/,,  <  /<  0.8/,  (band  1 ). 
0.8/,,  <  /<  1 2fp  (hand  2).  1.2/.  <  f<  2/  (hand  3).  and  2 f  <  /•  3/, 
(band  4).  The  KHH  n  2  run  corresponds  to  the  dissipation  pa¬ 
rameterization  detined  by  Rogers  el  al.  (2003)  and  based  on  KHH 
and  applied  in  ihe  WWATCH  code  For  the  peak  frequency  band, 
statistics  are  also  given  for  the  directional  spread  <rw 
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hourly  mean.  The  HRI)  winds  do  nol  cover  ihc  entire 
computational  domain  For  areas  falling  outside  the  do¬ 
main,  the  nearest  National  Data  Buoy  Center  (NDBC) 
wind  observation  is  used.  1  his  produces  some  non¬ 
physical  spatial  discontinuities  in  Ihe  wind  field,  but  these 


that  are  insufficient  to  resolve  small  storms  such  as  trop¬ 
ical  cyclones  (Tolman  and  Alves  2005).  Hurricane  waves 
do  share  many  similarities  with  more  usual  sea  slates 
(Young  2006),  but  the  high  winds  and  their  rapid  rotation 
are  particularly  challenging  for  numerical  wave  models,  ll 
is  thus  necessary  to  verify  that  the  new  source  functions 
perform  adequately  under  extreme  wind  conditions.  A 
simulation  of  Hurricane  Ivan  (Gulf  of  Mexico.  September 
2004)  is  chosen  for  this  purpose  because  it  was  extensively 
measured  (e.g.,  Wang  el  al.  2005)  and  hindeaslcd. 

Winds  for  this  simulation  are  based  on  3- hourly  grid- 
ded  surface  wind  analyses  created  by  National  Oceanic 
and  Atmospheric  Administration's  (NOAA)  Hurricane 
Research  Division  (HRD).  As  an  intermediate  step,  fields 
are  reprocessed  at  30-mm  intervals,  with  the  storm  po¬ 
sition  updated  at  each  interval  (thus,  scmi-Lagrangian 
inlet  polalion).  The  wind  speeds  are  reduced  by  factor  of 
1/1.1 1  to  convert  from  Ihe  maximum  sustained  gust  to  an 


0  8  I  12  14  1.6  18  2 

f/fp 


FUi.  13.  Model  data  comparison  it  NDBC  buoy  45007  (Lake 
Michigan),  (top)  Comparison  of  the  summed  spectral  density  vs  the 
frequency  for  the  duration  ot  the  simulation,  and  (bottom)  the 
correlation  coefficients  vs  Mu.  normalized  frequency  (see  text  for 
explanation). 
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arc  smoothed  in  the  wave  model  integration  and,  in  any 
event,  only  affect  weaker  wind  seas  far  from  the  storm 
center. 

Bathymetry  is  taken  from  the  Naval  Research  Labo¬ 
ratory’s  2'-resolulion  database  (DBDB2)  coarsened  to 
the  computational  grid  resolution  (0.13).  The  directional 
resolution  is  10°,  and  the  frequency  range  is  limited  to 
0.041 8-0.41 17  Hz.  The  model  was  applied  from  13  to  16 
September  2004.  Model  results  are  illustrated  in  Fig  14. 
The  models  were  validated  at  all  of  the  buoys  in  the  Gulf 
of  Mexico.  Results  from  buoy  42040,  where  the  waves 
where  largest,  are  shown  here. 

Model  runs  with  parameterizalions  BAJ,  TEST437, 
and  TEST441  give  very  similar  results:  close  to  the  ob¬ 
servations,  except  for  the  highest  waves  (//v  >  13  m  at 
buoy  42040),  where  TEST437  and  TEST44I  give  slightly 
smaller  values.  Results  with  the  TC  parameterization 
are  generally  lower  in  terms  of  fls  than  are  all  of  the 
parameterizalions  that  share  a  Janssen  type  of  input.  It 
appears  that  the  new  source  term  perform  similarly  to 
BAJ  and  is  able  to  reproduce  such  young  waves  and 
severe  sea  states. 

Because  the  wind  forcing  enters  the  wave  model  through 
a  wind  stress  parameterized  in  a  way  that  may  not  apply 
to  such  conditions  it  is  worthwhile  to  reexamine  some 
choices  made  above.  In  particular  the  surface  roughness 
w  as  allowed  to  exceed  0.002  in  TEST44I  b.  which  resulted 
in  better  estimates  of  lfs  l  ifting  this  constraint  shows 
that,  for  these  very  high  winds,  the  wind  sheltering  effect 
plays  a  similar  role  in  limiting  the  roughness  to  :1)m;ix, 
with  the  difference  that  it  tends  to  narrow  the  wind  input 
spectrum  (Fig.  7).  '1  his  narrower  wind  input  has  a  limited 
effect  of  the  wind  stress  and  //,,  but  is  has  a  noticeable 
effect  on  the  spectral  shape.  This  is  illustrated  by  the 
low-frequency  energy  that  appears  to  be  strongly  over¬ 
estimated  before  the  peak  of  the  storm  for  TFST441b. 
In  general  the  new  parameterizalions  provide  results 
that  are  as  reasonable  as  those  of  previous  parameleri 
zntions,  given  the  uncertainty  of  the  wind  forcing. 

5.  Conclusions 

A  set  of  parameterizalions  for  the  dissipation  source 
terms  of  the  wave  energy  balance  equation  have  been 
proposed,  based  on  known  properties  of  swell  dissipa¬ 
tion  and  wave  breaking  statistics.  ‘This  dissipation  in¬ 
cludes  an  explicit  nonlinear  swell  dissipation  and  a  wave 
breaking  parameterization  that  contains  a  cumulative 
term,  representing  the  dissipation  of  short  waves  by  lon¬ 
ger  breakers,  and  different  dissipation  rales  for  different 
directions.  These  dissipation  parameterizalions  have 
been  combined  with  a  modified  form  of  the  wind  input 
proposed  by  Janssen  ( 1991 ),  in  which  the  questionable 


IO  buoy 

•  T&O  (1996)  :  bias  = -0.80  m 

- BAJ  .  bias  =  -0.62  m 
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height  and  (bottom)  partial  wave  height  (for  f  <  0 .Oh  Hz),  during 
lhc  passage  of  Hurricane  Ivan,  at  buoy  42040.  The  10  m  discus 
buoy  was  capsized  by  waves  and  could  not  record  waves  after 
16  September. 

guslincss  parameter  r./r  has  been  reduced,  and  the  gen¬ 
eral  shape  of  the  wind  input  has  been  significantly  mod¬ 
ified.  The  resulting  source  term  balance  is  thus  markedly 
different  from  the  previous  proposed  forms,  with  a  near 
balance  for  very  old  seas  between  the  air  sea  friction 
term  that  dissipates  swell  and  the  nonlinear  energy  flux 
to  low  frequencies.  Also  the  wind  input  is  concentrated 
in  a  narrower  range  of  frequencies. 

For  younger  seas  the  wind  input  is  weaker  than  that 
given  by  Janssen  (1991 )  but  stronger  than  that  given  by 
Tolman  and  Chalikov  (1906).  However,  the  dissipation 
at  the  peak  is  generally  stronger  because  it  is  essentially 
based  on  a  local  steepness  and  these  dominant  waves  are 
(he  steepest  in  the  sea  slate.  As  a  result,  the  short  fetch 
growth  is  relatively  weaker  than  that  with  the  source 
term  combination  proposed  by  Bidlot  et  al.  (2007a)  (BAJ ). 
Flic  choice  of  parameters  tested  here  lends  to  produce 
broader  directional  spectra  than  those  observed  in  Lake 
Michigan  and  global  hindcnsls,  and  slanting  fetch  di¬ 
rections  that  are  too  oblique  relative  to  the  wind  (Fig.  2). 
In  this*  respect  (he  new  source  terms  arc  intermediate 
between  BAJ  and  TC, 

Another  defect  comes  from  the  definition  of  the  sat 
uration  level  used  to  define  the  breaking-induced  dissi 
palton.  Here,  as  in  the  work  by  van  der  Weslhuyscn  cl  al. 
(2007),  the  saturation  is  local  in  frequency  space,  whereas 
wave  breaking  is  naturally  expected  to  have  a  relatively 
broad  impact  due  to  its  localization  in  space  and  lime 
(Hasselmann  1974).  This  is  expected  to  produce  an 
overestimation  of  the  energy  just  below  the  peak,  and  an 
underestimation  at  the  peak  of  the  saturation  spectrum 


Sfptfmbkr  2010 


ARDHIJIN  ET  AI. 


1937 


TAHLI-.  Al.  Wind  wave  interaction  paranieiers  as  implemented  in  version 3.14  ofSMOM  in  ilicWW  A  lT'H  code,  and  values  used  in  ihc 
tesis  presenied  here.  In  WWATCH.  all  parameters  are  accessible  via  the  namelist  "SIN3,"  a  set  of  FORI  RAN  parameters  that  can  be 
defined  when  running  a  computer  code.  AH  of  these  parameters  are  included  in  version  3.14  of  WWATCH  Note  that  v„  is  a  switch  that,  if 
nonzero,  activates  the  calculation  of  5lMlt. 
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I  hose  effects  likely  contribute  to  the  persistent  over- 
estimation  of  low-frequency  energy  in  the  model. 

In  spite  of  these  defects,  the  new  parameterization 
produces  robust  results  and  clearly  outperforms  the 
Bidlot  et  al.  (2007a)  parameterization  in  global  hind 
casts,  whether  one  considers  dominant  wave  parameters. 
//*,  7 *m02'  and  lp  or  parameters  sensitive  to  the  high- 
frequency  content,  such  as  the  surface  Stokes  U ss  drift  or 
the  mean  square  slope.  At  global  scales,  errors  in  //v. 
and  Uss  are — on  average — reduced  bv  L5%,  25%,  and 
50%  relative  to  those  obtained  with  the  parameteriza¬ 
tion  by  Bidlot  et  al.  (2007a).  Global  and  regional  hind- 
casts  (North  and  Irish  Seas,  English  Channel,  Tuamotu, 
Lesser  Antilles,  and  more),  from  2(X)2  to  2010  at  least, 
are  available  online  for  further  analysis  (at  http://tinyurl. 
com/yetsofy). 

Another  important  aspect  of  this  study  was  the  vali¬ 
dation  at  regional  to  global  scales.  Tolman  and  Chalikov 
(1996)  included  verification  with  steady-state,  fetch-limited 
growth  curves.  Though  such  a  verification  is  a  useful  step, 
the  outcome  of  the  Lake  Michigan  hindeast  suggests 
that  such  verification  gives  no  indication  of  the  model's 
skill  in  real  subregional-scale  applications.  One  of  the 
parameterizations  proposed  here  (TEST441)  also  gives 
slightly  poorer  performance  for  young  seas,  which  is  not 
obvious  in  the  case  of  Lake  Michigan,  but  was  revealed 
by  hindcasts  of  Mediterranean  waves  (not  shown). 

Because  our  intention  was  only  to  demonstrate  the 
capability  of  new  dissipation  parameterizations  and  the 
resulting  source  term  balances,  we  have  not  fully  adjusted 
the  IS  parameters  that  define  the  deep  water  parame¬ 
terizations,  compared  to  the  roughly  9  used  by  Bidlot 
ct  al.  (2(K)7a)  The  results  presented  here  are  thus  pre¬ 
liminary  in  terms  of  model  performance,  which  is  why 
the  parameterizations  are  still  given  temporary  names 


like  TEST44L  As  illustrated  by  the  Hurricane  Ivan  hind- 
cast,  some  parameters,  such  as  arc  probably  un¬ 

necessary:  in  that  particular  case,  the  removal  of 
improved  the  results,  but  lor  global-scale  results  it  had  no 
impact  at  all  (not  shown). 

Because  five  of  the  extra  parameters  define  the  air  sea 
friction  term  that  produces  swell  dissipation,  and  two 
define  the  cumulative  breaking  term,  it  is  feasible  to  de¬ 
fine  a  systematic  adjustment  procedure  that  should  in¬ 
crease  the  model  accuracy  by  separately  adjusting  the 
swell,  wind  sea  peak,  and  high-frequency  properties.  In 
particular,  the  directional  distribution  may  be  improved 
by  making  the  dissipation  term  more  isotropic  (i.e., 
taking  <$,/  >  0.3)  or  modifying  the  definition  of  the  sat 
uration  parameter  B'  in  Eq.  (12).  In  a  forthcoming  paper 
we  shall  further  investigate  the  responses  of  the  wave 
field  to  varying  curients.  from  global  scales  to  regional 
tidal  currents.  From  that  research,  we  expect  that  wave 
steepening  will  produce  much  more  dissipation  due  to 
breaking,  as  envisaged  bv  Phillips  (1984). 

Obviously,  it  is  well  known  that  the  discrete  interaction 
approximation  used  here  to  compute  the  nonlinear  inter 
actions  is  the  source  of  large  errors,  and  further  calcula¬ 
tions  will  be  performed  using  a  more  accurate  estimation 
of  these  interactions  in  the  future. 
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El  Fluvialcs  (CETMF.F).  The  SHOM  buoy  deployments 
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Tablk  A2.  Dissipation  parameters  as  implemented  in  version  3.14  of  $HOM  in  the  WWA  ICH  code,  and  values  used  in  the  tests 
presented  here.  In  WWATCH,  all  parameters  are  accessible  via  the  namelist  "SDS.V  The  only  parameters  not  defined  in  the  present 
paper  are  Cu  and  Ch|,  which  act  like  switches  to  activate  the  B  AJ  parameterization  for  the  part  of  the  spectrum  with  saturation  below  and 
above  t  he  spectrum,  respectively  Most  of  these  parameters  are  also  included  in  version  3. 1 4.  except  for  Octl.  which  is  needed  for  TEST437 
and  TEST441  with  results  described  here  The  TEST405  can  be  run  with  version  3.14.  The  parameter  V/t,  is  a  switch  for  the  correction  or 
not  of  Br  into  B’r ,  when  Mu  1,  as  is  the  case  here,  the  correction  is  not  applied 


Parameter 

lEq- 

Variable  in  WWATCH 

WAM4 

BA! 

TEST405 

TEST437 

TEST44I 

FES  1443 

C\w 

(2) 

SDSCI 

4.5 

-2.1 

0.0 

— 

— 

r 

(2) 

WNMEANP 

-0.5 

0.5 

0.5 

— 

h  M 

(7) 

FXFM3 

2.5 

2.5 

2.5 

9  9 

A, 

(2) 

SDSDELTAI 

0.5 

0.4 

0.0 

— 

— 

(2) 

SDSDE1  TA2 

0.5 

0.6 

0.0 

— 

1  ds 

(13) 

SDS02 

0.0 

0.0 

-2.2  X  10  5 

— 

— 

Gr 

SDSLF 

1.0 

1.0 

0.0 

— 

— 

f  h. 

SDSHF 

1.0 

1.0 

0.0 

— 

(12) 

SDSDTH 

0.0 

0.0 

SO 

Br 

(12) 

SDSBR 

0.0 

n.o 

1.2  x  10  3 

9  x  10  1 

'Vu 

(IS) 

SDSBRFI 

0.0 

0.0 

0.0 

0.5 

2  x  Ccll 

(IS) 

SDSC3 

0.0 

0.0 

0.0 

2.0 

O.S 

•*/» 

(12) 

SDSCOS 

0.0 

0.0 

0.0 

2.0 

(13) 

SDSDC6 

0.0 

0.0 

1.0 

0.3 

0.3 

0.0 

A/i> 

SDSBMO 

0.0 

0.0 

1.0 

A/, 

SDSBM1 

0.0 

0.0 

0.2428 

m2 

SDSBM2 

0.0 

0.0 

1 3)995 

Mx 

SDSBM  3 

0.0 

0.0 

-2.5709 

— 

— 

M.i 

SDSBM4 

0.0 

0.0 

1.3286 

were  managed  by  David  German  with  invalauablc  help 
from  Guy  Amis.  Wave  data  were  kindly  provided  by  the 
U.S.  Navy.  ESA.  and  CNES.  and  the  many  in  situ  con¬ 
tributors  to  the  JCOMM  (WMO-IOC)  exchange  pro¬ 
gram  coordinated  by  Jean-Raymond  Bidlot. 

APPENDIX 

Parameter  Settings 

The  numerical  values  of  all  parameters  that  dclinc  the 
source  functions  arc  listed  in  Table  A1  for  the  wind-wave 
interaction  term  SHlm  and  Table  A2  for  the  wave-ocean 
interaction  term  3V*.  We  also  recall  that  the  nonlinear 
coupling  coefficient  (variable  Nl  PROP  in  WWATCH) 
is  set  to  2.78  x  107  in  all  cases,  except  for  the  two  pa- 
rameterizations  mostly  frequently  used  here,  with  Cn\  ~ 
2.5  x  1()7  in  TEST437  and  TEST441.  Although  the  best 
performance  for  most  parameters  is  obtained  with  the 
TEST441  settings,  its  underestimation  of  extreme  sea 
states  may  be  a  problem  in  some  applications  for  which 
TEST437  may  be  preferred.  A  full  tuning  of  the  model 
has  not  been  attempted  yet  and  it  is  possible  that  simple 
adjustments  to  /3max  C.cu,  rLU>  and  j„  may  produce  even 
better  results.  Finally,  these  parameters  have  been  mostly 
adjusted  for  deep  water  conditions  using  European  Centre 
for  Medium  Range  Weal  her  Forecasts  (ECMWF)  winds. 
Using  other  sources  of  winds  for  large-scale  applications 


may  require  a  retiming  of  the  wind  source  function,  which 
can  be  performed  bv  a  readjustment  of  0nv,.x. 
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